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Résumé 


Dans cette thèse, on étudie le nombre maximum de cycles limites qui bifurquent 
des orbites périodiques du centre linéaire & = y,y = —x perturbé par une classe 
généralisée d'équations différentielles de Liénard de la forme 


{ à = y — E(futr, y)y) — (fax, y)y), 
ÿ = —x — E(gu(x, y) + fa(x, y)y) — E(ga2(x, y) + far, y)y), 


Où fu, fai et gx sont des polynômes en x et y de degré {,n et m respectivement 
pour à = 1,2 et € suffisamment petit. 


On détermine ce nombre maximum pour cette classe des systèmes en utilisant la 
méthode de moyennisation du premier et du second ordre. De plus, Cette étude est 
illustrée par des applications. 


Mots clés : cycle limite, équation différentielle de Liénard, méthode de moyen- 
nisation. 


Abstract 


In this thesis, we study the maximum number of limit cycles for a class of poly- 
nomial differential systems that can bifurcate from the periodic orbits of the linear 
center & = y,ÿy = —x perturbed inside a class of the generalized Liénard polynomial 
differential systems of the form 


{ & = y — E(futr, y)y) — (far, y)y), 
ÿ = —x — E(ga(x, y) + fat, y)y) — (goal, y) + far, y)y), 


where f1;, fax and g»; are polynomials in + and y having degree /,n and m res- 
pectively for each à = 1,2 and € is a small parameter. 


The maximum number for this class of systems is determined, we use the avera- 
ging method of first and second order. This study is illustrated by applications. 


Keywords : Limit cycle, Liénard differential equation, averaging method. 
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Introduction générale 


En général par un système dynamique, on sous-entend tout système qui évolue 
au cours du temps. On représente cette évolution par des équations différentielles ou 
des applications. Ainsi, l'étude des systèmes dynamiques traite l’évolution temporelle 
des systèmes physiques, économiques, chimiques sans pour autant faire référence à 
la théorie sous jacente qui détermine leurs équations d'évolution. 

La notion des équations différentielles est aparut à la fin du 17°” siècle. A cette 
époque, les équations différentielles s’introduisent par le biais de problèmes d’origine 
mécanique ou géométrique comme par exemple : le mouvement du pendule circulaire, 
le problème du mouvement de deux corps s’attirant mutuellement suivant la loi de 
la gravitation newtonienne, … 


Un des principaux problèmes de la théorie des équations différentielles est l’étude 
des cycles limites : leur existence, leur nombre et leur stabilité. Un cycle limite pour 
un système différentiel plan, est une orbite périodique isolée dans l’ensemble de toutes 
les orbites périodiques. La notion des cycles limites stables est très importante car 
elle modélise des systèmes à oscillations auto-entretenues, voici quelques exemples : 
les battements du coeur, les vibrations d’un pont, ou des ailes d’un avion. Les cycles 
limites ont été introduits pour la première fois par H. Poincaré en 1881 dans son 
"Mémoire sur les courbes définies par une équation différentielle" [A0]. Poincaré s’est 
intéressé à l’étude qualitative des solutions des équations différentielles, c’est à dire 
aux points d'équilibre, aux cycles limites et leur stabilité, ce qui permet d’avoir une 
idée globale des autres orbites du système étudié. 


Un des théorèmes les plus importants de la dynamique non linéaire est le théorème 
de Poincaré-Bendixson qui affirme que dans une région bornée et compacte du plan, 
une trajectoire d’un système planaire converge vers un cycle limite ou un point 
critique. Pour la non-existence des solutions périodiques, il existe les critères de 
Bendixson et celui de Dulac, qui affirment sous certaines conditions que le système 
différentiel planaire n’admet aucune solutions périodiques. 
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En 1900, le mathématicien David Hilbert [T6] présenta, lors du deuxième congrés 
international des mathématiques , 23 problèmes "dont l’avenir attend la résolution 
grâce aux nouvelles méthodes qui seront découvertes dans le siècle qui commence". 
Le problème numéro 16 est de déterminer le nombre maximal et la position relative 
des cycles limites d’un système différentiel polynômial planaire de degré n. On note 
H, ce nombre maximal. En 1923, Dulac [12], proposa une démonstration prouvant 
que À, est fini pour tout n, mais sa démonstration comportait une erreur. La ré- 
solution de ce problème de Dulac a été faite de façon indépendante par Ilyashenko 
en 1991 et Ecalle, Martinet et Moussu en 1987 puis Ecalle 1992. Cette résolution 
permet de montrer que H, < c. Petrovsky, Landis en 1957, crurent trouver la va- 
leur de A; mais ils s’aperçurent d’une erreur dans leur propre démonstration (Landis 
et Petrovski, 1967) avant que celui-ci ne soit infirmé par un contre exemple de Shi 
en 1982, dans lequel un système quadratique a 4 cycles limites. Aïnsi, si H, est un 
nombre fini pour tout n, la seule chose que l’on sache est que H3 > 4 et H3 > 11 
(Jibin et Chunfu, 1985 ; Zoladek, 1995). Christopher et Lloyd ont donné en 1995 une 
borne inférieure au nombre H, : H, > n° log(n). 


En 1928, Liénard [24] un ingénieur français établit un théorème d’existence et 
d’unicité d’une solution périodique pour un système portant son nom de la forme 


nid 


ee. q 


où F(x) et G(x) sont deux polynômes de degré m et n respectivement. Notons que 
ce système inclue l'équation de Van Der Pol à + £(x° — 1)à + x — 0. Liénard 
montra que si G(x) = x et si F(x) est une fonction continue et impaire qui à une 
unique racine positive en x — a et qui est strictement croissante pour x > à, alors le 
système (1) possède un unique cycle limite. 


e En 1975, Rychkov [44] à prouvé que pour des polynômes F(x) impairs et de degré 
5 et si G(x) = x alors le système n’a au plus que deux cycles limites. 


e En 1977, Lins, De Melo et Pugh ont prouvé que si m = 3 et n = L,alors il 
n’y à au plus qu’un cycle limite. Ils ont de plus donné les conditions pour que 
ce cycle existe. Enfin ils ont conjecturé que si G(x) = x, il ne pouvait y avoir 
plus de E (#21) cycles limites (E(x) est la partie entière du réel x). Mais cette 
conjecture a été prouvée qu'elle est fausse pour n > 6. 


e En 1981, Zhang [50] a étudié le résultat de Rychkov (1975) à des fonctions non 
polynômiales : 
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— Si G(x) = x, si F(x) et sa dérivée f(x) sont continues. 


— De plus si la dérivée f(x) est pair et a deux racines a, a2 > 0 tel que a1 < à@2 
et Fa) > 0, F(a2) < 0. 


— Si f(x) strictement croissante pour x > @. 
Alors le système a au plus deux cycles limites. 


e En 1983, Xianwu a prouvé la conjecture pour certains cas où m = 4etn = 1. 
En 1990, Dumortier et Rousseau ont montré que H(3,1) = 1. 

En 1997, Duman et Li [3] ont prouvé que H(2,2) = 1. 

e En 2002, Wang et Jing [47 ont prouvé que H(3,2) = 3. 

e En 2012, Li.Ghengy et J.Llibre [8| ont prouvé que H(1,3) — 1. 


Cependent, certains chercheurs furent intéressés au nombre maximum de cycles 

limites noté A (m,n) appelé "Medium limit cycles" qui seraient crées par bifurcation 
d’un centre. Ainsi, il ont trouvé une borne inférieure de H(m, n). 
En effet, un aspect fondamental de l’analyse des systèmes dynamiques est la théorie 
de la bifurcation, introduite par Poincaré [A1] en 1892 pour décrire les changements 
qualitatifs des points d’équilibres d’une équation différentielle, obtenue selon une 
faible variation d’un paramètre, ensuite a été discuté par Andronov et Witt [3] et 
leurs collègues de travail commençant en 1930. 


En 2010, J.Llibre et al. [29] ont réussi d'obtenir des bornes inférieurs pour le 
nombre de cycles limites pour m,n > 1. Plus précisément, ils ont étudié le nombre 
maximum de cycles limites 27 (m,n) qui peuvent bifurquer des orbites périodiques du 
centre linéaire & = y,ÿy = —x perturbé par des équations différentielles de Liénard. 
D'autres travaux sont pointés sur la bifurcation des cycles limites pour des systèmes 
différentielles polynômiaux ayant un centre isochrone. (Voir chapitre 4). 


Dans cette thèse, on considère le problème de recherche des solutions périodiques 
pour un système différentiel dont le second membre dépend d’un petit paramètre €. 
Une réponse à ce problème peut être donnée en utilisant la méthode de moyenni- 
sation (averaging method), cette méthode est l’une des plus importantes méthodes 
de perturbations utilisée actuellement dans l’étude des cycles limites des systèmes 
dynamiques. Elle a été introduite en 1937 par Krylov et Bogoliubov et par Bo- 
goliubov et Mitropolskii [5] en 1961. Elle a été ensuite développée par Verhulst [46], 
Malkin (1956) et Roseau (1966) [43], en 2004 par Llibre et Buica [6]. D’autres formes 
et théorèmes de la méthode de moyennisation ont été démontrés ces dernières années 
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(voir [6|) et beaucoup d'articles ont été publiés concernant l'application de cette mé- 
thode. En utilisant la méthode de moyennisation, on réduit le problème difficile de 
la résolution des équations différentielles non autonomes à la recherche des racines 
positives d’une équation algébrique. 


En 2012, Llibre J. Valls.C ont étudié, en utilisant la méthode de moyenni- 
sation d’ordre un et deux un système d'équations différentielles de la forme 

& = y — e(gu(x) + fu(x)y) — E(gro(x) + fio(æ)y), 

ÿ = x — E(gn(x) + fa(x)y) — (gx) + fa(x)y), 


OÙ Qui, fai, Gui, Jui pour à = 1,2 de degré l,k,m et n respectivement et £ un petit 
paramètre, Ils ont montré qu’au plus 


À = max[y + [(m —1)/2],u + [1/2], Î(n — 1)/21+ [m/2], (k/21 + Im/2] — 1, 
[n—1)/2] +10 —1)/2]+1,1k/2] + [0 — 1)/2] 
(où u = min[n/2],[(k — 1)/2]) cycles limites bifurquent des orbites périodiques du 


centre linéaire & = y, ÿ = —x. 


En 2013, Llibre et Valls ont étudié le système différentiel polynômial plus 
général 


&=y— fi(x)y, 
ÿ = —x — ga(x) — fa(x)y, 


où f1(x) = Efu(x) + fio(x) + ES fis(x), go(x) = Ega(x) + gx) + ga(x), et 
fax) = Efa(x) + E foo(x) + fos(x) où fui, fui et gai de degré l,n et m respective- 
ment pour à = 1,2,3 et € suffisamment petit. 

Ils ont trouvé une borne supérieure du nombre maximum de cycles limites qui 
peuvent bifurquer des orbites périodiques d’un centre linéaire & = y, ÿ = —x, en 
utilisant la méthode de moyennisation d’ordre trois. 


Dans ce travail, on étudie les cycles limites des systèmes différentiels non auto- 
nomes en utilisant la méthode de moyennisation d'ordre un et deux. On s'intéresse 
à la recherche du nombre maximum de cycles limites qui bifurquent des orbites pé- 
riodiques du centre linéaire # = y, y = —x, perturbé par une classe généralisée 
d'équations différentielles de Liénard de la forme 


{ & = y —E(fulr, y}y) — (f(x, yy), 
ÿ = —x — E(gu(x, y) + fa(x, y)y) — (goal, y) + arf, y)y), 


OÙ fu, fa et ga; sont des polynômes en x et y de degré l,n et m respectivement pour 
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i = 1,2 et € suffisamment petit, ( voir chapitre 3). 
On note que ce système est plus général que celui étudié dans [30] et [31]. 


Cette thèse est subdivisée en quatre chapitres : 


Le premier chapitre est un rappel sur des notions générales. On introduira des 
définitions élémentaires tels que : le système dynamique, les points d'équilibre 
et leurs nature, la linéarisation au voisinage d’un point d'équilibre, les cycles 
limites, ainsi que des théorèmes sur l’existence et la non-existence des cycles 
limites. On introduira aussi un rappel sur la théorie de bifurcation. 


Dans le deuxième chapitre nous présentons la méthode de moyennisation du 
premier et second ordre, nous illustrons ceci à travers des exemples. 


Le troisième chapitre s'intéresse aux cycles limites des systèmes de Liénard, on 
détermine le nombre maximum de cycles limites qui bifurquent des orbites 
périodiques d’un centre linéaire perturbé par une certaine classe généralisée 
d'équations différentielles de Liénard et ceci en utilisant la méthode de moyen- 
nisation d’ordre un et deux. Cette étude à fait l’objet d’une publication inti- 
tulée : 

"S. Ellaggoune, S. Badi, On The upper bound of the number of limit cycles 
for a class of polynomial differential systems. Electronic Journal of Differential 
Equations. ISSN 072-6691. Volume 2016 (2016), Number 318, pp.1-22." 


Le quatrième et dernier chapitre est consacré à l’étude du nombre de cycles 
limites qui bifurquent des orbites périodiques d’un centre isochrone, nous in- 
troduisons la notion de ce dernier et celle d’intégrale première. 


— Chapitre 1 D 
Notions preliminaires 


1.1 Existence et unicité de la solution 


Soient U un ouvert de R x R” et f : U — R” une fonction continue. 


Définition 1.1.1. i) Une équation différentielle ordinaire (EDO) sur U est une re- 
lation de la forme : 


que l’on note brièvement 
& = f{t,x) où à = —. (Fi) 
ii) Pour (to,to) donné, un problème à valeur initiale associé à l’équation (1.1) 
est donné sous la forme 
{ à = f(tx), 


T(éo) = To: 


Définition 1.1.2. i) La fonction x(t) est dite solution de l'équation (1.1) sur un 
intervalle T CR si elle est définit et continûment dérivable sur T et si 


(t,x(t)) EU, VtelT et si x(t) satisfait la relation (1.1) sur J. 


(12) 


ii) Soit (to, to) € U donnée, la fonction x est dite solution du problème à valeur 
initiale (1.2) s'il existe un intervalle T1 contenant to tel que x est une solution 
de l’équation sur Ï et vérifie x(to) = %o. 

Définition 1.1.3. Considérons la fonction f(t,x) avec f : R"1 = R?, ft — to) < a, 

etre DCR? 

On dit que la fonction f(t,x) est Lipschitzienne par rapport à x, s’il existe K > 0 

telle que : 


[f(é, x) — fG, x2)| < K\x: = To], V(t,x1), (é, x) € [Éo — à, Lo + al x D. 


La constante K° est appellée constante de Lipschitz. 


1.2. STABILITÉ DE LA SOLUTION 


Théorème 1.1.1. On considère le système différentiel : 


r= f(x), TER" teR 


et on suppose que la fonction vectorielle f(t,x) est Lipschitzienne de rapport K° par 
rapport à x, uniformément en t € [-a, a]. Soit xo une donnée initiale, il existe une 
seule solution x(t) du système différentiel qui satisfait x(0) — xo et qui est définie 
sur l’intervalle [-—-c, c] avec c < min(a, +). 


Preuve. Cette solution satisfait l'équation intégrale 
t 
29 = 2(0)+ f Hu,r(u)du 
0 


on considère l’espace des fonctions continues y € C°([-a, a]) muni de la norme 
lui = marier-a,a y()Il. Soit L : C([-a, al) — C([-a, al) l'opérateur définit par 


LD = e(0) + fu vd 


cet opérateur satisfait 


et donc 
| LG) L L(y)| < K uv]. 


Si on pose c < min(a, +); on obtient que l’opérateur L est une contraction. Il possède 
donc un unique point fixe dans l’espace fonctionnel C([-—a, a]). Cet unique point fixe 
est une fonction continue qui est solution du système différentiel et ceci démontre 
l'existence et l’unicité cherchées. 


1.2 Stabilité de la solution 


La stabilité est l’un des aspects essentiels dans l’étude des systèmes dynamiques. 
Cette notion a été étudiée par Liapunov (1857-1918), voir le livre [39]. 


Définition 1.2.1. Soit le problème à valeur initiale 


SR (IT) 


Supposons que f satisfait les conditions du théorème d'existence et d’unicité de la 
solution. Une solution D(t) du système (1) telle que D(to) — Po est dite stable au 
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1.3. SYSTÈMES DYNAMIQUES 


sens de Lyapunov si : V € > 0, 3 0 > 0 telle que pour toute solution xt) de (1) 
dont la valeur initiale x(to) = ïo vérifie 


Ix(to) — Poll < 0 — [x(t) — P(t)| <e, Vt > to. 
Si de plus 
lim. led — &(01 = 0. 
alors la solution ®(t) est dite asymptotiquement stable. 


Remarque 1.2.1. Une solution ®(t) qui n’est pas stable est dite instable. 


Exemple 1.2.1. Soit le problème à valeur initial 


La solution de (1.3) est 
ét) =1. 


La solution de (+) telle que x(0) = xo est 
z(t) = (ro — 1)e* +1. 


On a : 
xt) — p(#)| = (to — 1)e | < |ro — 1}, VE > 0, 


il suffit de prendre Ô <E, 0 =Ee = @(t) = 1 est stable. 


1.3 Systèmes dynamiques 


Définition 1.3.1. Un système dynamique sur R” est une application 
U:R x R° — R” définie et continue telle que 


1) U(.,x): R — R” est continue. 


2) U(t,.): R° — R'" est continue. 
3) U(0,x) — 
4) U(t+s, U(t,U(s,x)), V,seR, Vre R. 


Proposition 1.3.1. Les systèmes dynamiques sont engendrés par des systèmes dif- 
férentiels. 


1.4. FLOT D'UNE ÉQUATION DIFFÉRENTIELLE 


Exemple 1.3.1. Soit le système différentiel linéaire 


DAT: 
ou (1.4) 


où À est une matrice constante, t € R* et x € R". La solution de (1.4) est donnée 
par 


a(t) = eÂxo 
Le système (1.4) engendre un système dynamique 
U:R'xR°—R", U(t,x) = eAx. 


En effet : 
1) On a Vt,T ER*,Vr e R': 


NUE + 7,2) — U(E,x)|| = flex — ex] < [ele x — ef] <eMAMeré — 7] II] 
< elAÏ (ITA 1)[x] — 0 quand r — 0, 


d'où U(.,x) est continue. 


2) On a Vt E R*,Vx,y Ee R': 
A A A 
AU, &) — UE, y) = lex — y = [lex — y < Mix — y, 
où M constante et puisque t est fixé d’où la continuité de U(t,.). 
3) On a U(0,x) = ex = Ir =. 


4) Vt,s eR*,Vr e R°, on a U(t+s,x) = eÂteAx = U(t,e#x) = U(t,U(s,x)). 
D'où le résultat. 


Définition 1.3.2. Un système dynamique U sur R° est linéaire si : 
U(E, ax + By) = à Ult,x) + 8 Ut,y), 
pour tout x,y ER", tEReta,BER. 


1.4 Flot d’une équation différentielle 


Soit le système non linéaire 


Bjr); MER. 


On note par /(x0) l'intervalle maximum d’existence de la solution @(#, x) du problème 


à valeur initiale 
ce f(x), xeR”, 


Z(0) = T5. de) 
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Définition 1.4.1. Soit E un sous ensemble ouvert dans R'" et f € C''(E), et soit 
pour x € E la solution D(t,xo) du problème (1.5) définie sur I(xo). Alors pour 
tE I(xo), l’ensemble des applications ®, définit par 


Pi(r0) = P(E, ro), 
est appelé le flot du système différentiel (1.5). 


d(xo) possède les propriétés suivantes : 


1. Oxo) est de classe C! 
2. Po(ro) = To 
3. Dtrs(ro) = D(d(r0)): 


Remarque 1.4.1. Le flot est dit autonome si f ne dépend pas explicitement du 
temps t, sinon ù est dit non autonome. 


Exemple 1.4.1. On considére le système : 


G= A, 
(0) = To 


QUEC 


La solution de ce problème est : O(t, xo) = exo avec 


= cost —sint 
| sint cost J’? 


donc le flot associé à ce problème est défini par : 


b : R? — R? 


To > (To) = exo = ( 


cost —sint 
sint cost La 


1.5 Théorie des systèmes différentiels non linéaires 
autonomes 


Dans cette section, on considère le système non linéaire autonome 


= f(x), æeR. (1.6) 
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AUTONOMES 


1.5.1 Points d’équilibre et linéarisation 


Définition 1.5.1. (Point critique). On appelle point critique (ou point d'équilibre 
ou point singulier ou point fixe) du système (1.6), tout point xo € R" qui vérifie : 


Remarque 1.5.1. Un point qui n’est pas critique est dit régulier. 


Définition 1.5.2. Un point critique x = a du système (1.6) est appelé attracteur 
positif s'il existe un voisinage V, de a tel que si : 


RE. li 0 =. 


—>+00 


Si un point critique x = a, a cette propriété quand t —> —, alors x = a est appelé 
attracteur négatif. 


Exemple 1.5.1. Soit le système 


T = —%, 
y — —2y, 
l’origine est le seul point critique pour ce système . 
On voit sur la figure 1.1 que le point critique (0,0) est un attracteur positif. 


— 


a 4 ie 6 — 6 — 6 — à — 


PS ll et ne à 


ee he — 5 — 2 gp en 


FIGURE 1.1 — Attracteur positif 
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Définition 1.5.3. (Linéarisation). Considérons le système non linéaire autonome 
(1.6), où xo est le point d'équilibre. 


Le système 


0j: . 
à = At, où A = Dj{ro) = 2 (to), 1<i,j <n, 
j 


(1.7) 
est appelé le système linéarisé de (1.6) en xo € R”. 
Exemple 1.5.2. Soit le système 
à = 22° — y, 
= (1.8) 
y=x+SYy, 


l’origine est le seul point critique pour ce système . 
La matrice jacobienne associée à (1.8) calculée en (0,0) est 


Dr00)= (1% ): 


T = —Y, 
Y = %, 


est le système linéarisé de (1.8) en (0,0). 


D'où le système 


Définition 1.5.4. Le point critique xo du système à — f(x) est dit hyperbolique si 
aucune des valeurs propres de la matrice jacobienne D f(x0) n’a de partie réelle nulle. 


Remarque 1.5.2. La linéarisation d’un système différentiel nous amène à l'étude 
de la nature des point critiques. 
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AUTONOMES 


1.5.2 Classification et nature des points d’équilibre 


Définition 1.5.5. On considère le système : 


Ge flo, ER, 


soit Xo son point critique. 


e Le point critique x) est appelé selle s'il est hyperbolique et si À = Df(xo) a au 
moins une valeur propre avec une partie réelle positive et au moins une valeur 
propre avec une partie réelle négative. 


e Le point x, est dit puits si toutes les valeurs propres de la matrice À = D f(x0) 
ont des parties réelles négatives. 


e Le point ro est une source si toutes les valeurs propres de la matrice À = D f(x0) 
ont des parties réelles positives. 


Exemple 1.5.3. Soit le système non linéaire autonome 
= 27-77, 
(1.9) 
= dr + 


Le système (1.9) a un seul point d'équilibre qui est l’origine (0,0), et le système 


linéarisé en ce point est 
z\ _/f20 TL 
ÿ) \OI y). 


Ce système a deux valeurs propres réelles de même signe À1 = 2 et À = 1. Alors le 
point critique (0,0) est une source. 


Pour un système plan à coefficients constants 


à = At, ze R’, 


où À une matrice carrée constante, l’origine (0,0) est le point critique. Considérons 
À et À2 les valeurs propres de la matrice À. On distingue les différents cas selon les 
valeurs propres Ài et À : 


(1) Si 1,2 Z 0 sont réelles et de signes différents, l’origine est une selle. Il est 
toujours instable. 
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| 
\ 
\ 
à 


FIGURE 1.2 -— selle 


(2) Si À; et À sont réelles de même signe, alors on à trois cas : 


— Si À1 < À2 < 0, l’origine est un noeud stable. 
— Si À > À2 > 0, l’origine est un noeud instable. 


— Si À = À2 = À, l’origine est un noeud propre, il est stable si À < 0 
et instable si À > 0. 


1.5. THÉORIE DES SYSTÈMES DIFFÉRENTIELS NON LINÉAIRES 


AUTONOMES 


ER un ne De De ne De 


4 0 Be 2e —Ù2ù 2e 0 te MN 


PR Re 20 pe 


[=] [=] 
prono eo 7 A4 ff 


Re ee er CT A 


PA 2 tt et mn à 


g 4 AS DT ET M M 


A a — 6 — 
6, 
oO (= | 


LS D ne 


Et Li ot out dur 1 : 
FARM SR 


— noeud stable 


FIGURE 1.3 


(3) Si À et À2 sont complexes conjuguées, alors l’origine est un foyer. Il est stable 


si Re(A12) < 0 et instable si Re(A12) > 0. 


PORT A AT I IT LE A Se ee 


A 4 
IN Re a à 


RSS 7 
RORS P  Ban e — —  r 


PALETTE Er er Er e 
LPS PTS EEE ee CU LR 
PIATNR  d'u dnts 

a | 

RE 

T1 

TIiN 

TX 

LEN 

RER 


6 LG eo a ne pe me où 


FIGURE 1.4 — foyer stable 
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AUTONOMES 


1.5. THÉORIE DES SYSTÈMES DIFFÉRENTIELS NON LINÉAIRES 


(4) Si À et À sont imaginaires pures, alors l’origine est un centre. Il est stable 


mais n’est pas asymptotiquement stable. 


EE D RE Tee un ue “45 “4e @à “à 


FOR A A IE LE AE 0 P —G Re ep ee nn 


FN IE AE SE LE Se D 0 of 66e #8 "0 "8 "8 


AN A LT NE AT ST Re ST 
AA A A AR I IT SE ro 


ge de ee es en 
En ee on en ee a 
PR en TS a Ve a 


BEBE OPSENT EF 
ee RRQ Be fe Les LR A 


“ d — nt nt A 7 71 
Pa Bo Les 


SR ST ET 
Se ni ot tt af a” de : 
ne 
9 A 


RS RS om — A A a AT mt met ne 


Centre stable 


FIGURE 1.5 — 


Exemple 1.5.4. Soit le système 


l’origine est le seul point critique pour ce système . 


—3. Alors le 


Ce système a deux valeurs propres imaginaires pures À1 = à et À 


point critique (0,0) est un centre stable. 
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1.5. THÉORIE DES SYSTÈMES DIFFÉRENTIELS NON LINÉAIRES 


AUTONOMES 
1.5.3 Plan et portrait du phase 
Définition 1.5.6. Soit le système planaire 
à = P(x,y), (110) 
ÿ = Q(x,y), 


où P et Q sont des polynômes en x et y à coefficients réels de degré n. 
Les solutions (xt), y(t)) du système (1.10) représentent dans le plan (x, y) des courbes 
appelés orbites ou trajectoires. 

Les points critiques de ce système sont des solutions constantes et la figure com- 
plète des orbites de ce système ainsi que ces points critiques représentés dans le plan 
(x, y) s'appelle portrait de phase et le plan (x, y) est appelé plan de phase. 


1.5.4 Théorème de Hartman-Grobman 


Le théorème de Hartman-Grobman (voir [39|) est un résultat important dans la 
théorie qualitative locale des systèmes différentiels. Il montre qu’au voisinage d’un 
point critique hyperbolique x, le système non linéaire 


à = f(x), (1.11) 


a la même structure qualitative du système linéarisé 
D Dis (1.12) 


Autrement dit, si l’origine est un point selle ou foyer ou noeud pour le système (1.12), 
alors le point critique x, sera respectivement selle ou foyer ou noeud pour le système 
(1.11). 

Cependant si l’origine est de type centre pour le système alors on ne peut rien 
dire sur la nature du point critique x, de (1.11). 


Exemple 1.5.5. Soit le système non linéaire 


à = x +27y, 
ÿ = —2y + 4x°y, 


donc l’origine (0,0) est le seul point critique de ce système, et le système linéarisé 


Df(0,0) = ( x “ } 
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en ce point est 


1.5. THÉORIE DES SYSTÈMES DIFFÉRENTIELS NON LINÉAIRES 
AUTONOMES 


Ce système a deux valeurs propres réelles de signes différents A1 = 1, A2 = —2. 
Alors le système linéarisé en (0,0) est un point selle instable. Donc le point critique 
du système non linéaire est aussi du type selle instable. 


Exemple 1.5.6. Soit le système non linéaire 


D=y=#, 
ÿ = —x — 3y°, 


donc l’origine (0,0) est le seul point critique de ce système, la matrice jacobienne en 


ce point est 
| 
pwo=(°, 1). 


Elle admet deux valeurs propres imagénaires pures À12 = +i. Alors le système linéa- 
risé on point critique (0,0) est un centre stable. Mais pour le système non linéaire 
on ne peut rien dire. 


1.5.5 Orbites périodiques et cycles limites 


Soit le système différentiel 


= JR) TeIRT. (1.13) 


Définition 1.5.7. (Orbite périodique). On appelle orbite périodique toute trajec- 
toire fermée o(t,x) du système(1.13) vérifiant 
O(t+T,x) = d(t,x). 
Le plus petit réel T > 0 qui vérifie cette égalité est appelé période. 
Proposition 1.5.1. Toute solution périodique entoure au moins un point d'équilibre. 


Définition 1.5.8. (Cycle limite). Pour un système plan, on appelle cycle limite 
une orbite périodique qui est isolée dans l’ensemble des orbites périodiques. 


Théorème 1.5.1. (Stabilité des cycles limites). C étant la trajectoire corres- 
pondante au cycle limite, toutes les trajectoires intérieures et extérieures voisines de 
C sont telle que : soit elles s’enroulent toutes en spirales autour de C pour t —> + 
où bien t —> —co. 


1. Le cycle limite est stable, si toutes les trajectoires intérieures et extérieures 
voisines sont attirées vers C'. 


2. Le cycle limite est instable, si toutes les trajectoires intérieures et extérieures 
voisines sont refoulées de C. 
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1.5. THÉORIE DES SYSTÈMES DIFFÉRENTIELS NON LINÉAIRES 
AUTONOMES 


3. Le cycle limite est semi-stable, si toutes les trajectoires intérieures et exté- 
rieures voisines sont attirées d’un coté est refoulées d’un coté. 


Exemple 1.5.7. Soit le système 


(1.14) 


E=2z-y—2r(x +y), 
ÿ = x + 2y — 2y(x° + y). 


En coordonnées polaires x — rcos(t), y = rsin(t) avec r > 0, le système (1.14) 
devient : 


F=2r(l=r), 
Ü= 1. 


d’où 
F=0—=r—-0 our = +l. 


Comme r > 0, on n'accepte que la racine positive r = 1. Donc, pour r = 1 on a la 
solution périodique (x(t), y(t)) = (cos(t + 60), sin(t + 6)), avec 0(0) = 6. 

Dans le plan de phase il y a un seul cycle limite dont l'équation x? + y? = 1 et 
d'amplitude r = 1. Voir figure (1.6) 


FIGURE 1.6 — Cycle limite du système (1.14). 


Remarque 1.5.3. Les cycles limites apparaissent seulement dans les systèmes dif- 
férentiels non linéaires. 
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1.6. EXISTENCE ET NON-EXISTENCE DES CYCLES LIMITES 


1.6 Existence et non-existence des cycles limites 


Théorème 1.6.1. {Poincaré-Bendixson). Soit le système plan suivant 


. = {(,9), (1.15) 


y = g(x, y), 


supposons que f,g sont des fonctions de classe C? sur E, où E est un sous 
ensemble ouvert de R?, le système a une orbite 7 telle que l'orbite positive 
7+(p) = {®(p,t),t > 0} passant par le point p est contenue dans un sous ensemble 
compact À de E. Alors on est dans l’un des trois cas suivants : 


e Soit y,(p) tend vers un point d'équilibre. 


P 


-(p) tend vers une orbite périodique. 
e Soit y,(p) est une orbite périodique. 


e Soit y: 


Une conséquence très importante de ce théorème est quand À ne contient pas de 
points critiques où dans ce cas, on est certain de l’existence d’une solution périodique. 


La difficulté dans l’application de ce résultat est bien la construction de la région 
(elle correspond souvent à un anneau) À. Pour cela cherchons deux cercles (C1) et 
(C2) telle que (Ci) soit à l’intérieur de (C2). 

La recherche des cercles (C1) et (C2) ressemble à la recherche de la fonction de 
Lyapunov v(x,y) = 2? +y?, 24+y?, xt + y. 

Soit (C) : u{x,y) =? +y = c 


a >0:(C):r+y =Q 
CG >0:(C2):7°+7y = ©, telle que «à < © 


Soit (C) : u(x,y) =x + =c>0. 
Si pour tout point p € (C') 


On a 


a M 
MP dy” 


— Toute trajectoire qui passe par le point p se dirige vers son extérieur. 
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Exemple 1.6.1. Montrons que le système 


=y-2 +5, 
ÿ=-x—Y +, 


admet au moins une solution périodique dans une région de R?. 


On a le point (0,0) est le seul point d'équilibre, la matrice jacobienne en (0,0) 


est : 
1 1 
a-(i) 


cette matrice admet deux valeurs propres complexes À = 1 +, À2 = 1 —i, donc 
(0,0) est un foyer instable pour le système linéarisé = il reste foyer instable pour le 
système non linéaire d’aprés le théorème de Hartman-Grobman. 


Cherchons l’anneau À : 
soit v(x,y)= 2? +y =0c 


dv dv dv 
HE a? ty — x + x) + 2y(—x — y° + y) 


= 2[(2° + y) — (2° + y]. 


On a pour 


CEE no 
2 + y < x? + y 


d 
donc si on prend (Ci) : x? + y? = « avec 0 < & < 1 et T > 0 


— toute trajectoire coupant (C1) se dirige vers l'extérieur. 
On a pour 


Br + y” > 2 
2(2° + y") > (x? + y)? > 2(x° + y?) 


dv 
P >2 —<0 
OUT Co dt 


— toute trajectoire coupant (C2) se dirige vers l’intérieur. 

Alors d'aprés le théorème de Poincaré Bendixon il existe au moins une solution 
périodique dans l’anneau À formé des deux cercle (C1) et (C2) puisque cet anneau 
ne contient pas de points critiques. (Voir la figure) 
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os oi ne nf nf af a dd 2 
Ca us nt et Ai A 7 A 2 
Sn nt nt a a dt 


mn ue ee me au ae at tt nn 


1 « 
FOR UT NO IT LR EE ne ee ne D D 
RIT IT DE RE AE EE Re 0 8 ee a a 
IT UT AT AT IT ET A ee qe pe ee nn eg 
ee eh — 2 a A ee 2 


PIE A I TE RE LT mm 
ST eV Te ne ne TS 

SP a De, Po a, ne 
Re Be a Se si 6 4 


Théorème 1.6.2. (Critère de Bendixson). Soit le système plan 


a 


ÿ = g(x, y), 


et soit F = (f,g)! € CE) où E est une région simplement connere dans R?. Si 
la divergence du champ de vecteur F (notée VF) est non identiquement nulle et ne 
change pas de signe dans E, alors ce système n’a aucune orbite fermée entièrement 
contenue dans E. 


Preuve. 


Supposons que (1.15) possède une solution périodique de période T noté 7, elle 
correspond à une courbe (orbite) fermée dans E. 
D'’aprés la formule de Green on a : 


Î/ (OF, y) F M ardy = | Jay = ue 
G=intrieur de T Ôx 0y Tr 


OT ECO 


= | LAGE(E), yO))g GE), y) — gx), y) f(x), y(t))]dt 
0. 
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Ceci est une contradiction, puisque divF est non identiquement nulle et ne change 
pas de signe, c’est-à-dire, soit négatif soit positif, donc l’intégrale ne peut pas être 
nulle. Alors il ne peut pas avoir de solution périodique contenue dans Æ. 


Exemple 1.6.2. Considérons le système suivant 
à = 2xy — 2y", 
ÿ=2 —ÿ — 2. 
Soit F = (2xy — 2y{,x? — y? — gs). On calcule la divergence du champ de vecteur 


F', on obtient 


OP — _ 0 D TS 
divF = VF = 2 (22y A) + y x°y) 


y y y = Sr y e 0. 


D'où, d’après le critère de Bendixson ce système n’a aucun cycle limite dans R?. 


Théorème 1.6.3. (Critère de Dulac). Soit F = (f,g)! € C\(E) où E est une 
région simplement connexe dans R?. S'il existe une fonction B € C\(E) telle que 
V(BF) est non identiquement nulle et ne change pas de signe dans E, alors le 


système 
à = f(x, y), 
ÿ = g(x,y), 


n'admet aucune orbite fermée entièrement contenue dans E. 


Preuve. 


Supposons qu'il existe une solution périodique 7 de période T contenue dans £. Elle 
est représentée par une trajectoire fermée dans le plan. 
Soit G l’intérieur de cette trajectoire d’après la formule de Green on a : 


} L CL | 1) re 1: 7 


: B{f(n@ — g(x,y) dé 


1.6. EXISTENCE ET NON-EXISTENCE DES CYCLES LIMITES 


Ceci est une contradiction, puisque divBF est non identiquement nulle et ne change 
pas de signe, c’est-à-dire, soit négatif soit positif, donc l’intégrale ne peut pas être 
nulle. Alors il ne peut avoir de solution périodique contenue dans Æ. 


Exemple 1.6.3. Soit le système suivant 


T'Y, 
ÿ = —ax — by + ax? + By°. 


Soit B(x, y) = be 28 avec F = (y, ax — by + ax? + By?)”. Calculant la divergence 
de BF : 


g 6) 
PB) 2 Ja ie (= a) 
— —28be y Fe be 28 Se 26be y _ _ pe 28% ES 


Donc, il n'existe pas de solutions périodiques pour ce système. 


e Méthode pour obtenir les fonctions de Dulac 
La liste des fonctions B données ci-dessus ne sert qu’a un nombre limité de systèmes. 
On discute quelques situations plus générales. 


On considère le système 


É = fi(x, 2) (1.16) 


To = fo(x1, T2) 


1. Premièrement, proposons une autre fonction positive ou négative K(x1,22) qui 
s’annule seulement sur un ensemble de mesure nulle telle qu’elle satisfait 


(AB) | O(RB) 
Ox: | 0x) 


—= K(x1, To) 


2. Maintenant choisissons À pour obtenir quelques résultats. 
Alors prenons K(x1,%2) = c(x1,72) X B(x1,22), avec c est une fonction positive ou 
négative qui s’annule uniquement sur un ensemble de mesure nulle, et substituons 
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cette relation dans l’équation précédente pour obtenir l’expression suivante : 


0B Of. Of: 
has ‘i n ie + — — c(x1, T2) X B(x1,2%2) (LT) 
On peut réecrire cette . 
Of. Of2 
D + ge = B(c(a2) _ Lu. F —) (1e 


3. Trouvons la solution de cette équation aux dérivées partielles quasilinéaire du pre- 
mier ordre. 

4. Finalement, voyons si B est une fonction de Dulac pour le système (1.16) jus- 
qu'ici, on a construit une méthode qui nous permet d’écarter les orbites périodiques 
associées au système d'équations différentielles dans le plan telle qu’elle est résumée 
dans ce qui suit. 


Théorème 1.6.4. Pour le système d’équation différentielle (1.16), une solution B 
du système associées (pour quelques fonctions c qui ne changent pas de signe 
et s’annulent seulement sur un sous-ensemble de mesure nulle), est une fonction de 
Dulac pour le système dans toute région À simplement connexe contenue dans 


D\{B7"(0)}. 


Corollaire 1.6.1. Pour le système d'équations différentielles (1.16), si (pour 
quelques fonctions € qui ne changent pas de signe et s’annulent seulement sur un 
sous-ensemble de mesure nulle) a une solution B sur D telle que B ne change pas 
de signe et s’annule sur un sous-ensemble de mesure nulle, alors B est une fonction 


de Dulac pour le système sur D. 


Remarque 1.6.1. Avec le choix de K dans l'étape 2, on obtient une équation linéaire 
(1.16) qui est simple à manipuler. 


Remarque 1.6.2. Avec le choix de c est utilisé pour simplifier l'équation (1.17), on 
peut toujours prendre c une constante non nulle. 
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1.6. EXISTENCE ET NON-EXISTENCE DES CYCLES LIMITES 


e Applications 


Exemple 1.6.4. Considérons le système 


d1 = tix(1 — cos to) 
Lo —= 32Ÿ%2 


d’aprés l'équation (1.18), on obtient : 


0B 0B 
zix(1 — cos to) — + 3220 — = Blc(x1, 2) — (2xix(1 — cos x) + 3x°)] 


0%: Ox) 
qui est une équation différentielle aux dérivées partielles quasilinéaire. 


En appliquant la méthode des caractéristiques, on obtient le système associé suivant 


dr, = xix5(1 — cos xo)dt 
dis = 3x? rodt 


dB = Blc(x1,x2) — (2r1x2(1 — cos x2) + 3x°)]dt 
des deux premières équations, on élimine le paramètre t on obtient 
3x Ÿtodti = rir5(1 — cos r2)dto. 
Résolvant cette équation, on obtient la première caractéristique 


1 
3T1 + Lo Sin To + COS Lo — 5% —;C;. 


Prenons c = 31° dans la dernière équation et on multiple l’équation par x1, 
on trouve 


r1dB = —B(2x%x5(1 — cos x2))dt. 
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Exemple 1.6.5. Considérons le système suivant 


T1 — —2%1%9 
2,2 2 2 
La = 2122 COS 21 + (1 + x )To + 22 


Utilisons l'équation (1.18), et prenons c = 1 + x, on obtient 


0B 0B 
—2tit2— + (xŸx5 cos ri + (1 + x? )xe + 25) — = —B(2rŸr2 cos ti). 


OT: Ôx 


Supposons que B dépend uniquement de x1 on a 


Tètitr = — B(2x°x%2 cos ti), 
i 


alors la solution est 


= exp( / T1 COS T1dT1). 
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1.7 Théorie de bifurcation 


Les systèmes d'équations différentielles paramétrées peuvent avoir différents com- 
portements asymptotiques (tendre vers un équilibre, un cycle limite, .…) en fonction 
des valeurs de leurs paramètres. Ils existent certaines valeurs pour lesquelles le com- 
portement du système passe d’un état qualitatif à un autre. Ce changement d'état 
qualitatif est une bifurcation et la valeur associée du paramètre est appelée valeur 
de bifurcation. 


La théorie des bifurcations des champs de vecteurs a pour but de décrire les 
modifications du portrait de phase des champs de vecteurs qui dépendent différen- 
tiablement d’un paramétre y 


ë= f(x n), (1:19) 
oùuER, fEC'(E), E un ouvert de R”. 


Considérons l’exemple suivant : 
soit le système différentiel paramétré 


_—_—. 
L A (1.20) 


Y = —x + y. 


Le système (1.20) a un seul point d’équilibre (0,0). la matrice associée est : 


1 y 
elle admet les valeurs propres À12 = y +1. 


e Si < 0, alors l’origine est un foyer stable. 

e Siu > 0, alors l’origine est un foyer instable. 

e Si = 0, cette matrice admet une pair de valeurs propres imaginaires pures, donc 
l’origine est un centre. 


Quand ce paramètre change de signe (c-à-d passant du négatif au positif), le 


système change de stabilité et passe de stable à instable. 
On dit, qu'il y a une bifurcation et dans ce cas y = 0 est la valeur de bifurcation. 
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1.7. THÉORIE DE BIFURCATION 


m>0 


Il existe deux types de bifurcations : les bifurcations locales et les bifurcations 
globales. Les bifurcations "locales" sont appelées ainsi car elles peuvent toujours 
être identifiées lors d’une linéarisation du système au voisinage de la solution. Le 
critère de détection utilisé dans le cas des bifurcations locales concerne les valeurs 
propres du Jacobien qui intervient au premier ordre dans la linéarisation. Par contre, 
les bifurcations "globales" ne font donc pas forcément intervenir le voisinage de la 
solution, ici, les linéarisations locales autour de la solution ne seront donc d’aucune 
aide. C’est pour cela que ces bifurcations sont appelées globales. Ci dessous, un 
exemple d’une bifurcation locale. 


Exemple 1.7.1 (Bifurcation selle-noeud). Soit le système différentiel suivant : 


: 2 

T=X + 

| * (1.21) 
ÿ = —y 


Selon le signe de 4, on distingue trois cas : 


e Sir < 0 dans ce cas le système admet deux points d'équilibre de coordonnées : 


(V0) et (1,0). 


La matrice jacobienne du système s'écrit 
J U 


J(x,y) = ( + ) 
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Au point (—/—u,0), la matrice jacobienne s'écrit 
2/=ù 0 


Cette matrice admet deux valeurs propres réelles négatives À1 = —1 et ]9 = —2\/—y, 
d’où le point d'équilibre (—,/—u,0) est un noeud stable. 


Au point (4/—u,0), la matrice jacobienne s'écrit 
2/— 0 
J( 2 , 0) — ( 0 # _] ) 


Cette matrice admet deux valeurs propres réelles et de signes opposés À; = —T et 
À2 = 24/—u, d’où le point d'équilibre est un point selle instable. 

Donc, quand  < 0 le système admet deux points d’équilibres, le premier du type 
noeud et le deuxième du type selle. 


e Si = 0, le système (1.19) devient 


i = 7? 
Y = —y 


L'unique point critique du système est l’origine. La matrice jacobienne associée s'écrit 


J(0,0) = ( . L ) 


Donc le point d'équilibre (0,0) est non hyperbolique. L'étude de la première équa- 
tion montre que x = 0 est un point d'équilibre non hyperbolique et la deuxième 
équation montre que y = 0 est un point d'équilibre asymptotiquement stable. 


e Si > 0, dans ce cas le système (1.21) n'admet aucun point d'équilibre car 
trek. 
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Cette bifurcation est une bifurcation locale appelée "selle-noeud" correspond à 
l’apparition de deux points d'équilibre, l’un stable (un noeud) et l’autre instable{un 
point selle) et au point de bifurcation u = 0 les deux points d'équilibre disparaissent et 
un point non hyperbolique apparait, qui en u > 0 lui même disparait : cette bifurcation 
correspond donc au changement de nombre de points d'équilibre et leur stabilité. 


Dans cette thèse, on s'intéresse à un autre type de bifurcation locale appelée 
la bifurcation d’un centre linéaire. On considère un système différentiel plan 
paramétré (qui dépend de €), pour lequel il admet un centre à l’origine quand € — 0, 
on étudie s’il y à naissance d’orbites périodiques lorsqu'on varie le paramètre €. 
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| Chapitre 2 ; 
Méthode de moyennisation 


La méthode de moyennisation est l’une des plus importantes méthodes de pertur- 
bations utilisée actuellement dans l’étude des cycles limites des systèmes dynamiques, 
cette méthode donne une relation quantitative entre les solutions d’un système diffé- 
rentiel périodique non autonome et celle de son système différentiel moyenné lequel 
est autonome. 


Dans ce chapitre, nous présentons la méthode de moyennisation d’ordre un 


et deux appliquée aux équations différentielles contenant un petit paramètre € de la 
forme 


Défaite) 


oùteR,xeR. 


Définition 2.0.1. On considère la fonction f : R x R° —> KR” qui est continue 
et T-périodique en la première variable t. On appelle fonction moyennée f° de f la 
fonction définie par 


fa = 2 j fe æ)dt. 


Définition 2.0.2. Soit f : R x R° —> R” une fonction continue. On dit que f a 
une moyenne notée f° si la limite 


existe et 


= fr, z)dr — f(x) < ko(T), V(t,x) ER XxXR" 


2.1. THÉORÈME DE MOYENNISATION 


avec k > 0 et o : [0,+oo) — [0,+oc) est une fonction strictement décroissante, 
continue et bornée telle que o(T) — 0 quand T — oo. 
2.1 Théorème de moyennisation 
Considérons le problème à valeur initial suivant 
 =ef(t,x) +eg(t,x,e), x(0) = eo. (2.1) 


On suppose que f(t,x) est T-périodique en t et on introduit la moyenne 


LT 
P= x [fu dt 
0 
Considérons maintenant le problème 


ÿ = ef" (y), y(0) = vo. (2.2) 


Théorème 2.1.1. Soient les problèmes aux valeurs initiales (2.1) et (2.2) avec 
xz,y,t%o E DCR”, t > 0. Supposons que 


1. f,g et D,f sont continues et bornées par une constante M indépendante de € 
dans 
[0,+co) x D. 


2. g est Lipschitzienne en x € D. 
3. f(t,x) est T-périodique ent, T est indépendante de €. 
A. y(t) est contenue dans un intérieur de D. 
Alors, on a x(t) — y(t) est de l’ordre O(E) pendant un temps d'échelle 1/€. 


Où la solution y(t) représente une approximation de x(#). 


Lemme 2.1.1. (Gronwall). Supposons que pour to <t <to+T 


o(t) < AU — to) + à | o(s)ds + 03, 


où @(t) est une fonction continue, @(t) > 0 Vto <t <to+T, 01,02 et 03 sont des 
constantes avec 01 > 0, 0 > 0 et 03 > 0. Alors 
02 


te) < (2 + sgpeñt-to © 
1 


di 
pour to <t <to+T. 
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Preuve du Théorème (2.1.1). 


Les hypothèses (1) et (2) nous assurent l’existence et l’unicité de la solution des 
problèmes aux valeurs initiales (2.1) et (2.2) sur l'échelle du temps £. 
Posons : 


(6,9) = L (Fe,2) = Fu))ds. 


to 


On a 
lu(é, y) <2MT pour t > to et y E D. 


Maintenant, nous introduisons 
2(t) = y) +eu(t, y(t)). 
Comme y(t) € D, Vt > to, on a l'estimation 
Ie(E) — JO < x) — 20 +120 — YO 


< xt) — 2461 +ellu(, yO)| 
< |x(t) — z(t)|]| + 2e MT. 


Notons que 
t'dx dz 
x(t) — z(t) — se 
0-20 = [QD 
Calculons 
dr dr. - dy 
am dr) +e tte) dt 


CO CO) 


d 
Remplaçons a par ef°(y), on aura 


æ - - On 
R = g(t,x(t),€) — TE y) Fu) — eff y()) — Ef(E (0). 
On a 


If (y) < M et ECO < 2MT. 


Grâçe à la continuité lipschitzienne de f nous avons 


LAC, 2()) — FC y) < LI20) — y) 
< eLu(s, y(4))| 
< 2LMT. 
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Donc il existe une constante k telle que 
LR] < ke? 
C’est clair que 
dé. dé 
20 < | | - <a 
le 2 < [1 - Gé 


<e | f(x) — FC, 2) + kellds 


to 
t 
= «| LF(E, 2 (6) — FE, 2(6))IIds + KeŸ(t — to) 
to 
t 
EL | rte) — 2(0)lds + ke — 1) 
to 
D'’aprés le lemme de Gronwall 
k k 
Ix(é) — z(t)|]| < ere L(t — to) — ET: 
Par conséquent 
fs k 
Ir — y < LE — to) — 2 + 2MT) 
Si eL(t—to) est borné par une constante indépendante de €, on aura l’approximation 


x(t) = y(t) — O(E) quand € — 0. 


31 


2.2. THÉORÈME DE MOYENNISATION DU PREMIER ORDRE 


2.2 Théorème de moyennisation du premier ordre 


On considère le système différentiel 


i=eF(t,x) += G(t,x,e), x(0) = x, (2:3) 


où x € D CR”, D est un domaine borné et t > 0. Supposons que F(t,x) et G(t,x,€) 
sont des fonctions T-périodiques en t. 
Le système moyenné associé au système (2.3) est 


g =Ef{y), y(0) = %o, (2.4) 


où 


PQ) 7 | Fe ds. 25) 


Théorème 2.2.1. Soit le système (2.3), on suppose que F,G, D,F, D?F and D,G 
sont continues et bornées par une constante M dans [O, 00) x D avec —Eo < € < €o. 
Supposons aussi que F et G sont T-périodiques en t, où T indépendante de £. Alors 
on à : 


(1) Si le point p est un point critique pour le système moyennisé (2.4) tel que 
det (D; f°(p)) # 0. (2.6) 


Alors pour |e| > 0, suffisamment petit, il existe une solution T-périodique x(t, €) 
du système (2.3) telle que x(0,€) —> p quand € — 0. 


2) Si le point critique y = p du système moyenné (2.4) est hyperbolique, alors 
P que y P y 
pour |e| > 0 suffisamment petit, la solution périodique correspondante x(t,€) 
du système (2.3) est unique, hyperbolique et de même stabilité que p. 


Exemple 2.2.1. Soit L'équation de Van Der Pol 


ÿ+z=e(l-x)à. (2.7) 


L'équation (2.7) peut s’écrire sous la forme d’un système à deux équations différen- 
tielles 


Ÿ =, 

y — —x+e(1—x), 
en coordonnées polaires (r,0) où x = rcos(0), y = rsin(4) avec r > 0 , ce système 
devient de la forme 
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ñ = er sin*(0) (1 — r” cos”(8)) , 2.8) 
Ô = —1+ecos(0)sin(#) (1— 7° cos”(8)). | 
Le système est équivalent à 
dr —_— _ p2 2 “2 2 
= € (1 — r° cos” (8)) sin°(8) + O(e?). (2.9) 


On note que l'équation (2.9) est sous la forme standard (2.3), ainsi on peut appliquer 
la méthode de moyennisation en prenant 


z=7r, t=0, T=2»r et F(r,0) = —r (1-7 cos(8)*) sin*(8). 
On calcule l’équation (2.5) on obtient 


27 
Poe m2 r (1— 7° cos”(4)) sin? 940 = L (r? —4). 
27 Jo 8 
f°(r) a une unique racine positive r = 2. Puisque (df°/dr) (2) = 1, par l'hypothèse 
(1) du théorème l'équation de Van Der Pol (2.7) a pour |e| Z 0 une orbite 
périodique d'amplitude r = 2. De plus, puisque (df°/dr) (2) = 1 > 0 par l'hypothèse 
(2) du théorème (2.2.1), ce cycle limite est instable. Voir la figure (2.1). 


pe a a “a a NX 


= 
> 
Ts 
D 
= 
R" 
à 
\ 
\ 


PAPA A 


FIGURE 2.1 - Cycle limite instable de l'équation (2.7) pour € = 0.01. 
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Exemple 2.2.2. Soit le système 
T = —y+e(-x+xy), 
ÿ=z+e(-y+yz+2), (2.10) 
è=e(y+x +uyz—1), 


en coordonnées cylindriques x = rcos(0), y = rsin(0) et z = z, le système (2.10) 
devient 


d 

T = € (r° cos®(0) sin(0) + sin(0)2? — r +rz —rzcos°(6)), 

dô 1 

1 1 + . (rz cos(8) sin(9) + cos(0)2° — cos”(4)r° + r° cos" (@)), 
_. = € (rsin(8) +r° cos (0) + rzsin(0) — 1). 


Pour écrire le système (2.10) sous la forme standard pour appliquer la méthode de 
moyennisation on considère 0 comme nouvelle variable indépendante. D'où 


_ = € (r° cos?(@) sin(@) + sin()2? —r+rz-—7rz cos! (8)) | 

(2.11) 
d 
D =# (r sin(0) + r°? cos?(0) + rzsin(0) — 1) | 


Le système (2.11) est de la forme 
dr 


LE 2 
dô eFi(r, b, z) 7 O(E j; 


dz 9 
d0 — EFPo(r, 0,2) + O(E h 


où F1, F2 sont périodiques de période 27, Calculons maintenant le système moyenné, 
on obtient 


27 
fn = 3 [ Ftn62)d8= 372), 
27 Jo 9 
1 [2 1, (2.12) 
fa(r, z) — x | P(r,0,z)dû —= 97 — I. 


On résout le système (2.12), on obtient les deux racines 
(r1, 2) — (V2, 2) ; (T2, 22) = (-v3 2) | 


Comme r > 0, la solution (V2, 2) est la seule qui fournisse un cycle limite. On 
vérifie maintenant que le déterminant calculé en cette racine est non nul. 


of ôf nn 
_ Se D _ =g=T] 2 __ 
D(r1, 21) = det ôf, BR, det 2 . 2 | = = —1 £0. 
Or Oz 
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D'après le théorème (2.2.1), le système (2.10) a pour || suffisamment petit un 
seul cycle limite. 
Ofi Of 


De plus la matrice FA ÿ, calculée au point (Va 2) admet les deux valeurs 
Ôr  z 


HO ‘0x 
propres +1 (réelles de signe différent). Donc le point singulier est un point selle qui 
est toujours instable. Alors, le cycle limite est instable. Voir la figure (2.2). 


FIGURE 2.2 - Cycle limite instable de l’équation (2.10) pour € = 0.001. 


2.3 Théorème de moyennisation du deuxième ordre 


Lorsque le champ de vecteurs est suffisamment différentiable, on peut passer à la 
méthode de moyennisation d'ordre deux. 
considérons les deux problèmes aux valeurs initiales 


i=ellt,z)+e Pt,x) + Ps(t,x,e), æ(0) = x, (2.13) 


où F1, F2 : [0,+oo) x D —> R”, F3 : [0,+co) x D x [0,€0] —> R” sont des fonctions 
continues, T périodiques par rapport à t et D est un ouvert de R”. 


ÿ = ef (y) + [F (y) +9 ()], y(0) = To, (2.14) 
où f0, f1% et g° sont des fonctions moyennées correspondantes à F, f, et F; respec- 
tivement. 


39 


2.3. THÉORÈME DE MOYENNISATION DU DEUXIÈME ORDRE 


Théorème 2.3.1. Soit 


où 
t 
HET) = | [F(s, ) — f(x) ds + z(x), 
0 
et z(x) est une fonction de classe C! dont la moyenne est nulle. 


Supposons que 
OF 


(a) ——,P2 et F3 sont Lipschitziennes relativement à x et continues dans leur do- 
maine de définition. 
(b) |Bi(t,x,€)| est bornée par une constante M positive dans [0, %) x D x (0,€0]. 
(c) Test indépendante de €. 
1 
(d) y(t) € D pendant un temps d'échelle —. 
E 
Alors 
a(t) = y(t) + eyi(t, y(E)) + O(E?), 
1 
sur l’échelle du temps —. 
€ 
Corollaire 2.3.1. Si les hypothèses du théorème (2.3.1) sont satisfaites et de plus 
f°(y) = 0. 
Alors 
1. Sip est un point d'équilibre du système moyenné (2.14) tel que 


5 (Fu) + (y). 40! (215) 


Y—P 


Alors, il existe une solution T-périodique @(t,£) du système (2.13) telle que 
@(t,E) — p quand e —> 0. 


2. Si (2.15) est négative, la solution périodique o(t,E) du système (2.13) est 
asymptotiquement stable pour n suffisamment petit. Si (2.15) est positive, cette 
solution est instable. 


Exemple 2.3.1. Considérons le système 


a | 


2.16 
g=z+e (y°x + y — 2y + 6y° — y°). ) 


En coordonnés polaires (r,0) où x = rcos(@), y = rsin(#) avec r > 0, on obtient 
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= e(—4r* cos*(0) sin(@) + r° cos*(0)) + &?(r* cos(0)sin(8) — r° cos*(0)sin(0) 
+ 2r cos?(8) — r$ cos’ (8) + r°cos$(@) — 2r + 6r° + r$ cos(8) — r° — 3r$cos®() 
+ 3r° cos/(0) — 12r° cos” (0) + 3r° 0 — 3r° cos*(0) + 6r° cos“(@)), 
Ô = 1+e(4r? cos?(8) — 4 cos{(0)r? — sin(8)r cos?(0)) + &2(— cos(0)r{ sin(0) 
+ 2sin(8) cos*(0)r* — 6r° cos*(0) 7 ) + sin(@)r° cos*(0) + sin(0)r° cos°(0) 
— 2sin(0r° cos{(0) — cost (8)r? — sin(8) cos” (0)r* + r? cos?(0) + 2sin(8) cos*(0)r* 
+ 6r° cos(0) sin(0)). 


Considérons maintenant 0 comme la variable indépendante, on obtient 


_. = eFi(r,0) + eP(r, 0) + O(E*), (2.17) 


où 
F\(r, 0) = —4r° cos*(0) sin(0) + r? cos*(8), 
P(r,0) = —16r° cos’(8) sin(0) + 16r° cos”(0) sin(0) + r° cos” (0) sin(9) — r° cos*(@) sin(0) 
+ r% cos(8) sin(@) + r° sin?(0) cos°(0) — 8r“ cos®(0) sin?(0) — r° sin?(0) cos“(8) 
— 25 sin( 0) cos®(@) + 2r° sin?(4) cos?(0) — 6r° sin?(0) cos?(0) + r6 sin?(0) cos(6) 
— r° sin?(0) — 2r sin?(0) + 6r° sin°(0). 
En appliquant maintenant le théorème (2.3.1), on calcule la fonction moyennée 
1 27 
PA) = 37 | F(r.0)d8 = 


On peut passer à la méthode de moyennisation d'ordre deux, on calcule 


g°(r) = ) “ht 0)d0 


27 
et à 
= fo. / 
Fri 2x J, —— (7,0): | Fi(r,s)ds| dû 
où 
A (r. 0) = —8r cos*(8) sin(0) + 2r cos*(8), 
et 


/ { 2 
| Fr, s)ds = —r* + r° cos“ (0) + 37" cos" (8) sin(0) + à sin(0). 
0 
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Par conséquent, 


res [co [res +rt0)] à 


5 4, 90 
=r|-— —r —1). 
r ( 16” + ru ) 
Pour trouver les cycles limites, on résout l’équation (2.18) et on obtient deux ra- 
cines positives r1 = 54/(90 + 1061), r2 = 54/(90 — 1061). 


Etant 


(2.18) 


nous obtenons 

CAT 0 

—(f (71) + g°(r1)) = —26.25844945 £ 0, 
pa 

et 


EF) + g°(r2)) = 1.85844941 Z 0. 


Donc il existe deux cycles limites, et d'après l'hypothèse (2) du corollaire (2.3.1) 
on à 


d d 

Un) + g°(r1)) = —26.25844945 < 0 et TU (r2) + g°(r2)) = 1.85844941 > 0. 
id Tr 

Alors, le premier cycle limite est stable d'amplitude r1 — ë (90 + 1061) et le 


deuxième cycle limite est instable d'amplitude r2 = à (90 — 1061) pour € suffi- 
samment petit. Voir la figure (2.3). 
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/ 
PENSE LISE D. 
$Ÿ 00 MN SES À A 
Dot, D DMMAA A A A 


a of af ar er a eh CR KO 


4 Ok x du % = es 2 
BRL LL = 
& 1 EN D Re deb 


À nm % 5 = ps SR 


FIGURE 2.3 — Deux cycles limites du système (2 pour € = 0.001. 
Exemple 2.3.2. Considérons le système 


{ = -y+e (x —4r°y) + (y — xy), 


2.19 
g=x+e (ufr + 9x — 4y + 6y° — y°) . ie) 


En coordonnés polaires (r,0) où x = rcos(4), y = rsin(0) avec r > 0, on obtient 


F = e(—4r* cos*(8) sin(0) + r? cos? (8)) + e°(—r? cos” (8) sin(0) — r° cos” (8) 
+ 4r cos? (0) + r? cos(@) — 4r + 4r°cos*(0) + 6r? + r$ cos(0) — r° — 3r$cos*(0) 
— 5r° cos“(0) — 12r° cos” (0) + 3r° cos°(8) — r$ cos’/(0) + 2r° cos° (0) + 6r° cos“(8)), 
Ô = 1+e(4r? cos (8) — 4 cos{(0)r? — sin(4)r cos?(0)) + £2(—2 cos*(0)r° sin(0) 
— r{sin(0) cos(8) + r° cos°(0) sin(@) — 2 sin(0)r* cos” (0) + sin(0}r 
— 6sin(@)r? cos*(0) + 3sin(0) cos*(@)r* — r sin(0) — 4sin(4) cos(8) + r cos(@) + sin(0) cos?(0)r° 
— r cos”(0) + 6r° cos(@)sin(0)). 


cos” (0) 


Considérons maintenant 0 comme la variable indépendante, on obtient 


D 2 <P(r,0) + 2 Fa(r,0) + O(e) (2.20) 


où 
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F\(r, 0) = —4r* cos*(0) sin(4) + r° cos*(8), 

Fs(r, 0) = —16r° cos’ (0) sin(0) + 16r° cos” (8) sin(8) + r? cos”(0) sin(@) — r? cos?(0) sin(0) 
+ 76 cos°(0) sin?(0) — 8r*sin?(0) cos” (0) — 2r° cos! (0) sin?(9) — 27° sin?(@) cos*(0) 
+ 37° sin?(0) cos*(9) — 62r° sin?(8) cos?(0) + r? sin?(0) cos(0) + r° sin?(4) cos(0) 
— 4r sin?(9) + 6r° sin?(0) — r° sin?(8). 


En appliquant maintenant le théorème (2.3.1), on calcule la fonction moyennée 


27 
PA) = 37 | F(r.0)d8 = 


On peut passer à la méthode de moyennisation d'ordre deux, on calcule 


1 27 
0 _— — 
g (r) = x [ Pi(r, 0)d0 
et , : 
f'(r) = x) EG 6) , F\(r, sy dô 
A Jo à 0 
où 
F 
A (r. 0) = —12r° cos” (0) sin(0) + 2r cos*(0), 
et 


É 1 2 
| Fr, s)ds = —r* + r° cos“ (8) + 37” cos"(8) sin(0) + Se sin(0). 
0 


Par conséquent, 


Fr) + g(r) = sf E (r, 0) - Î[ Fi(r, s)ds + Fi(r, 0) | dô 
=r (-rt+ Tr _ ) 


Pour trouver les cycles limites, on résout l'équation (2.21) et on obtient deux ra- 
cines positives T1 = 1,72 = 24/2. 
Etant 


(2.21) 


d 10 0 _ 9 4 27 2 
Ar) + gr) = = grt + Dr 2, 
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nous obtenons 


LPO) + g)(r1)) = 3500000000 £ 0, 


et 


LPO) + gra) = 28 7 0. 


Donc il existe deux cycles limites, et d’après l'hypothèse (2) du corollaire (2.3.1) 
on à 


d d 

, (fr) + g°(r1)) = 3500000000 > 0 et a UV" (r) + g°(r2))) = —28 < 0. Alors, 
Tr r 

le premier cycle limite est instable d'amplitude r1 = 1 et le deuxième cycle limite est 

stable d'amplitude r> = 2V2 pour € suffisamment petit. Voir la figure (2.4). 


CAS 7 A de de de HU 
EF # ET AT ET AT EE, 6, 0 OÙ À 
ÉÉE # LS #7 EME 
“ FE P a À SON OÙ RO OX 
ÉÉÉ PA a 4e 5%. 9 Où RO NX À 
VV (a NONNRRA 
É É É V # » » V6 R Rk 
té É LA ER 
111811: \liii 
Tiitiis LE 
t + + Y à 
VV à \ k à A 1 1 1 
RRIIR D | 
RARES A 4 À 
RER RTS 4 A # f 
Yu fs A A À 
Vu a a a a na DH A + 4 A 
UNS USA RD bb 7 p ft À A A 
w D 'E." 2 del + #7 #7 VAT A 


FIGURE 2.4 — Deux cycles limites du système (2.19) pour € = 0.001. 
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2.4 Un autre théorème de moyennisation du pre- 
mier ordre 


Théorème 2.4.1. On considère le problème de bifurcation des solutions T-périodiques 
du système différentiel de la forme : 


ë(t) = Pot,z) +eFi(t,x) +e*P(f,x,e), 52) 


où € € (—E€0,€0) pour £0 suffisamment petit. Les fonctions F6, F : R x Q — R” 
et PB :RXx x (—E60,€0) — R” sont des fonctions de classe C?, T-périodique 
ent et Q est un ouvert de R”. Supposons que le système non perturbé 


&(t) = Fo, x), (2.23) 


a une sous variété des solutions périodiques de dimension Kk. 
Soit x(t,z) la solution du système non perturbé telle que x(0, z) = z. 
La linéarisation du système non perturbé le long de la solution périodique 
x(t,z) s'écrit 
ÿ = D, Pot, ft, z))y. (2.24) 


Notons par M,(t) la matrice fondamentale du système différentiel linéaire (2.24). 
Supposons qu'il existe un ensemble ouvert V avec C'L(V) € Q, tel que pour chaque 
z € CL{(V), x(t,z,0) est T-périodique, où x(t,z,0) est la solution du système non 
perturbé avec x(0,2,0) = 2. L'ensemble C'L(V) est isochrone pour le système 
(2.22), c’est à dire il est un ensemble formé seulement par des orbites périodiques, 
toutes ayant la même période. 

On note par £ : R° x IR? — RÉ Ja projection de R" sur ses k premières coor- 
données, c’est à dire 


A Et HE Un 


et par Et: REXRT É — IR F Ja projection de R' sur ses n—k dernières coordonnées, 
c’est à dire 


Et (1, me) = nm 


Alors, on a les résultats suivants : 


Théorème 2.4.2. [36) Soit V un ensemble ouvert et borné de RF, et soit 
8: CL(V) — R'À une fonction de classe C?, supposons que : 
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(i) Z = {2 = (a, B(a)),a e CL(V)} CA et pour chaque za € Z la solution x(t, za) 


de (2.23) est T-périodique. 


(ii) Pour chaque za € Z, il existe une matrice fondamentale M.,(t) de (2.24) telle 
que la matrice MLY(0) — MZY(T) a dans le coin supérieur droit une matrice 
k x (n—k) nulle, et dans le coin inférieur droit la matrice A,(n—k) x (n—k) 


avec det(As) £ 0. 
On considère la fonction F : CL(V) — R* 


F(a) = ( | : ML(t)Fi(t, xt, za))dt). (2.25) 


S'il existe a € V avec F(a) = 0 et det((dF/da)(a)) £ 0, alors il existe une solution 
T-périodique p(t,E) du système (2.22) telle que 4(0,€) — 2 quand € — 0. 


Nous donnons maintenant un résultat quand n = k. 


Théorème 2.4.3. [36) Soit V un ensemble ouvert et borné avec CL(V) € tel que 
pour chaque za € C'L(V) et pour chaque z4 € Z, la solution x(t, z4) est T-périodique. 
Considérons la fonction F : C'L(V) — R” 


FN 7 | MH E(6, dé za))dt. (2.26) 


S'il existe a € V avec F(a) = 0 et det((dF/da)(a)) £ 0, alors il existe une solution 
T-périodique p(t,E) du système (2.22) telle que @(0,e) — 2 quand € — 0. 


Théorème 2.4.4. [36] Supposons que n = 2m. Soit V un ensemble ouvert et borné 
de R" et soit B: CL(V) — R" une fonction de classe C?, supposons que : 


(i) Z = {za = (a, B(a)),a e CL(V)} CA et pour chaque z4 € Z la solution x(t, za) 


de (2.23) est T-périodique. 


(ii) Pour chaque z, € Z, il existe une matrice fondamentale M,,.(t) de (2.24) telle 
que la matrice M5'(0) — MI'(T) a dans le coin supérieur droit la matrice 
Afm x m) avec det(A,) Æ 0, et dans le coin inférieur droit une matrice 
m x m nulle. 
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On considère la fonction F : C'L(V) — R” 
1 [T 
F()=Et (> 1 MDP (Er za))dt). (2.27) 
TJ ” 


S'il existe a € V avec Fa) = 0 et det((dF/da)(a)) £ 0, alors il existe une solution 
T'-périodique p(t,E) du système (2.22) telle que 4(0,€) — 2, quand € — 0. 


Exemple 2.4.1. On considère l'équation suivante 
E —&+-—x=Ee(2+sin(t})(x? +4x°). (2.28) 


On écrit l'équation différentielle du troisième ordre (2.28) comme le système diffé- 
rentiel du premier ordre suivant 


T —=Y 
ÿ=2 (2.29) 
2=x—y+2+e(2+sin(t))(x? +47). 


L'origine est l’unique point singulier du système (2.29) lorsque £ = 0. La partie li- 
néaire du système (2.29) avec £ = 0 à l’origine est 


D: LL ‘0 
A= | 0 O0 !1 
1 12 1 


Les valeurs propres de la matrice sont +i et 1. 
On va faire un changement de variables linéaire 


(x, y, 2)7 = B(X,Y, ZYT, 


telle que dans les nouvelles variables (X,Y, Z,V), le système (2.29) avec € = 0 a sa 
partie linéaire égale à sa forme normale de Jordan, c.à.d. (X,Y,Z)7 = J(X,Y,Z)T, 
la forme normale réelle de Jordan de la matrice A est 

—1 


0 0 
J= T1 0 0 
0 O0 I 


On a 
BAB !=J= BA-JB=0, 
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d’où 
1 —1 0 
B — O0 —1 1! 
1 O0 1 


Par la transformation linéaire inversible (X,Y, Z)T = B(x,y,2)7, c.à.d. 


pe x X à 
Y sl vis TT SRB vi. 
Z 2 _ É 
on trouve 
X =i— Ÿ, 
Y = -ÿ+3, (2.30) 
Z=à+2. 
On a 
T X 
le] € |: 
z Z 
on obtient 
x X—-Y +7 
y |=-=| —-X-Y +Z |, (2:51) 
Z —X+Y +7 
on remplace (2.29) et (2.31) dans (2.30), on trouve 
X=-Y 
Y=X+Er(X,Y,Z,t) (2:82) 
Pret 
où 


F = F(X,Y,2,t) = F(x,y,2,t). 
Pour € = 0, la solution du système (2.32) est 


X(t) X9 cos(t) — Yisin(t) 
Y(t) | = {| Ycos(t) + Xosin(t) 
Z(t) Zoe 


On utilise la notion introduite dans le théorème 2.4.1, et d'aprés le système (2.32), 
on à 
r=(X,Y,2), F(x,t) = (—Y,X,2), Fi(x,t) = (0, F,F) et Py(x,t,e) = (0,0, 0). 
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Soit xt, Xo, Yo, Zo,€) la solution du système (2.32) telle que 
æ(0, Xo; Yo; Zo; €) = (Xo; Yo, Zo). 


IL est clair que, le système non perturbé (2.32) avec £ = 0 admet un centre à l’ori- 
gine dans le plan (X, Y). Les solutions périodique de ce centre sont x(t, X0, Yo, 0,0) — 
(X(t),Y(t), Z(t)) telle que 


X(t) X9 cos(t) — Yosin(t) 
Y(t) | = | Yocos(t) + Xosin(t) 
Z(t) 0 


Notons que toutes ces orbites sont périodiques de période 27. 


1. Pour notre système, V et a du théorème 2.4.2 sont 
V = {(X,Y,0),0 < X? + Y? < p} pour certains p arbitraires 
et a = (Xo, Yo) € V. 


2. La matrice fondamentale M(t) du système non perturbé (2.32) est 


cos(t) —sin(t) 0 
Mt) = | sin(t) cos(t) 0 
0 0 e! 
D'autre part, un calcul simple donne 
0 0 0 
M71(0) — M (27) = | O0 0 0 | 
0 O0 1—-e 


d'où 1e ?7T 20. 


Nous avons montré que toutes les hypothèses du théorème 2.4.2 sont vérifiées. 


Par conséquent, nous allons étudier les zéros a = (X0, Yo) € V des deux premiers 
composantes de la fonction F(a) donnée par 


F(a) — ex J MOAtaft, za))dt), É(t1, To, T3) = (x1, To), 


c. à. d. 


F(a) = (Fi(a), F(a)), 
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d’où 
= 2 f" sin(t)F(x(t, Xo, Yo,0,0),t)dt, (2.33) 
_ ip. | RECE Y(t ) _X( x6+7 on AO 
= f” cos(t)F(x(t, Xo, Yo, 0, 0), t)dt, (2.34) 
- up. cos(t)F( 0 w : -À€ “are TE —2,pat 


On pose F(a) = 


On intégre (2.33) et (2.34), on obtient 


l 3 3 3 3 
A) 30% = 370 Xo _ 100 = et = 30 


1 
FX, Yo) = 5 (A0 — M = Xi} 


Si Fi(X0, Yo) = F2(X 0, Yo) = 0, on trouve 
1 1 
Ne 
(XVe) = (@) 


On à HF, F) 
1,92 : 
pt lents) _ ne F0 


Alors, pour € € [—E€0,€0] avec €o > 0  . petit, il y a une solution isolée 
2x-périodique x(t,E) de l’équation différentielle (2 telle que 
1 
x(0,e) — 0, #(0,€) — 3 &(0,€) — 0, 


quand € — 0. 


AT 


- Chapitre3 
ycles limites pour une classe 


généralisée des systèmes de 
Liénard 


Dans ce chapitre, on étudie les cycles limites qui bifurquent des orbites pério- 
diques du centre linéaire & = y,ÿy = —x perturbé par une classe généralisée d’équa- 
tions différentielles de Liénard. 

Cette étude à fait l’objet d’un article publié dans le journal " Electronic Journal of 
Differential Equations" : 


S. Ellaggoune, S. Badi, Upper bounds for the number of limit cycles of polyno- 
mial differential Systems. Electronic Journal of Differential Equations. ISSN 072-6691 
Volume 2016 (2016), Number 318, pp. 1-22. 
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3.1 Historique 


Nous allons étudier dans ce chapitre des équations différentielles du deuxième 
ordre. d 
y =EF(y,y) où y=yia)et =. 
Ces équations peuvent se transformer en un système de 2 équations différentielles du 
premier ordre. 
L'espace des phases est alors un plan. La forme générale d’un tel système est : 


(3.1) 


Le premier modèle physique publié dans la littérature qui, transformé en un système 
du type (3.2), admette un cycle limite est l'équation : 


ae) 1 += ” 


établie par RAYLEIGH (1945) et qui modélise les oscillations d’une corde de violon. 


e Dans les années vingt, Balthasar van der Pol, un ingénieur hollandais, étudiait 
les propriétés électriques des tubes à néon (Van der Pol, 1922). A cette époque là, les 
oscilloscopes n’existant pas encore, il surveillait l’évolution de son circuit en écoutant 
les changements de tonalité dans un combiné téléphonique. Il modélisa les charges 
et décharges du tube par l’équation qui porte maintenant son nom : 


2 


d°x dx 


2 
— — 1 = 0, 3.4 

de +e(x Le sn” (3.4) 

a | sn dy 

On voit que si on dérive l’équation (3.3) par rapport à t et que si l’on note x — . 


on retrouve l’équation (3.4. Ces deux équations sont donc équivalentes. 

Quelques années plus tard Van der Pol (1927) étudiait le même circuit électrique 
mais en régime sinusoïdal forcé. Lorsqu'il variait la fréquence du courant, il entendait 
dans son combiné la tonalité changer par à-coups (le circuit se stabilisait donc sur 
la fréquence externe). Mais de temps à autre il remarquait quelque chose d’étrange, 
un comportement inexplicablement irrégulier : " on entend souvent au téléphone un 
bruit irrégulier avant que la fréquence ne saute à la valeur immédiatement inférieure" 
(Van der Pol, 1927). Son tube à vide devait vraisemblablement traverser une période 
de chaos transitoire avant de se synchroniser sur la fréquence externe. 

Plus tard, en Angleterre, Mary Lucy Cartwright et John E. Littlewood poursui- 
vront les travaux de B. Van der Pol sur les oscillateurs forcés. 
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e Liénard (1928), un ingénieur français, établit un théorème d'existence et d’uni- 
cité d’une solution périodique pour une classe générale d'équations dont fait partie 
l’équation : 

dx dx 


de - fa) kæz=0 (3.5) 


Liénard a transformé l’équation (3.5) en un système de deux équations du premier 


x 
ordre en posant — = 2 : 


dt 
dx dz 


nr en 0 (3.6) 


En fait dans l’énoncé de son théorème, Liénard utilise un autre système, équivalent 


à (3.6), obtenu en posant z = y — F(x), où F(x) = f} f(r)dr : 


dx dy 


x ?— ); FINE (3-7) 


Le plan (x, y) est alors appelé le plan de Liénard. L’équation (3.5) est appelée l’équa- 
tion de Liénard et on se réfèrera au système comme au système de Liénard. 

Le théorème de Liénard établi que si F(x) est une fonction continue et impaire 
qui à une unique racine positive en x —= a et qui est strictement croissante pour 
x > a, alors le système possède un unique cycle limite. 


e Levinson et Smith (1942), ont suggéré de généraliser le système (3.7) : 


dx dy 
Fa = U— (x) ; dt — —G(x), (3.8) 
qui se met sous la forme : 
d?x dx 


et qui est connu sous le nom de système de Liénard généralisé. 


e Le problème fondamental lié au système est le nombre de solutions pé- 
riodiques isolées (i.e. cycles limites) qui peuvent exister simultanément. Imaginons 
que le système décrit le mouvement d’un oscillateur. Si ce système a un cycle 
limite globalement attracteur alors l’oscillateur va évoluer, aprés un régime transi- 
toire, selon un mouvement périodique. Le point important est que la période de ce 
mouvement sera la même quelle que soit la condition initiale. 

On voit donc que la présence d’un cycle limite peut être une propriété importante 
d’un système surtout si ce système est une horloge. 
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3.2 Perturbation du centre linéaire 


Nous allons considérer la perturbation indépendante du temps où la fonction 
perturbée est le centre linéaire ( _ } 


0 


(i)-(2)-(): 


qui est bien le système de Liénard. 

Pour € = 0, il n’y à pas de cycles limites, mais seulement des trajectoires périodiques 
du type centre. 

La perturbation introduite transforme le centre du système linéaire en un foyer et 
transforme occasionnellement certaines trajectoires périodiques du centre en cycles 


Si on choisit f(x) — ( 


) le système étudié sera : 


limites, La méthode de moyennisation permet de connaître le nombre de cycles limites 
qui apparaissent. 


3.3 Bifurcation des cycles limites des orbites pé- 
riodiques d’un centre linéaire perturbé par une 
classe généralisée d’équation différentielle de Lié- 
nard 

Nous allons étudier le nombre maximum de cycles limites qui peuvent bifurquer 


des orbites périodiques du centre linéaire & = y, y = —x, perturbé par une classe 
généralisée d'équations différentielles de Liénard de la forme 


{ = y — fix, y)y, 


RE 0 


avec f1(x, y) = Efu(r, y) +e*fialx, y), fax, y) = Efn(x, y) +E* faa(x, y) et gox, y) = 
Ega(t,y) + E*gro(x, y), Où fu, fai et gx des polynômes en x et y de degré l,n et m 
respectivement pour à = 1,2, et € suffisament petit. On note que ce système est plus 
général que celui étudié dans [31]. 


Nous considérons le système (3.11) 


ÿ = —x — E(gn(x,y) + fa(x,y)v), 
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Où f11: 921 et fai des polynômes en x et y de degré {,m et n respectivement, et € 
suffisament petit. 


Notre résultat est le suivant : 


Théorème 3.3.1. Pour € suffisamment petit, le nombre maximum de cycles limites 
pour la classe généralisée du système de Liénard (3.11) qui bifurquent des orbites 
périodiques du centre linéaire à = y, ÿ = —x est 2] utilisant la méthode de 


moyennisation d'ordre un. 


Ensuite, nous considérons le système 


& = y — E(fu(x, y)y) — E(fiolx, vu), (3.12) 
ÿ = 2 — E(gn(r, y) + fn(x,y)y) — E(ga(x, y) + f(x, y)y), 

où f11 et f12 de degré l ; gx1 et g22 de degré m, et f21, f22 de degré n et € suffisamment 
petit. 

On a : 


Théorème 3.3.2. Pour € suffisamment petit, le nombre maximum de cycles limites 
pour la classe généralisée du système de Liénard qui bifurquent des orbites pé- 
riodiques du centre linéaire à = y, ÿ = —x en utilisant la méthode de moyennisation 
d'ordre deux est : 

À = max{ hi, À + 1, 3 +2}, 


E(m) + O(1) 1 
2 } | 2 


E(m) + O(n) — 1 
2 1 


A2 = max{O(n)—1, | 9 [AI 9 IA 9 LE 
Br) = 2,260) + 00 3, jOU) + ED 8, 
— EC un — =) 


e Oi) est le plus grand nombre impair inférieur ou égal à 1. 
e Eli) est le plus grand nombre pair inférieur ou égal à i. 
e |.] désigne la partie entière. 
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Preuve du Théorème (3.3.1). 


On utilise la méthode de moyennisation d’ordre un. Pour cela posons 


l 


Ju(x, y) nu r aijit Y?, fa(x, y) —. aijat y ; gai( T 4) s bijoz y ; 


i+j—=0 i+j=0 i+j=0 


En coordonnées polaires (r,0) où x = rcos0, y = rsin0, r > 0, le système (3.11) 
devient 


n 


= —e( Ÿ_ aijar 911 cos’ (9) sin}*?(0 +Y bij27"T1 cos’(0) sin?*1(0) 


i+j=0 i+j=0 


l 
+ à. Guar  cos (0) sin/*1(8)), 
i+j=0 


ô=-1- 2e ( . aij2r TT cos +1(0) sini*1(0) + >. bij2r T1 cos *1(0) sin/(0) 


i+j=0 i+j=0 


l 
— >. aijar 171 cos’ (0) sin/+2(9)) |. 
i+j=0 


(3.13) 


Considérons maintenant Ÿ comme nouvelle variable indépendante, le système (3.13) 
s'écrit sous la forme standard du théorème de moyennisation d’ordre un 


d 
= EF(r,0) + O(e?), (3.14) 
d0 
où 
F(r,6) = 2 aij27"7)+1 cos! (0) sin/*?(0 4e bar cos! (0) sin/+!(0) 
i+j=0 i1j=0 
l 
Fe D ajari titi cosi+! (0) sin/*1(0). 

i+j=0 

Maintenant, on calcule la fonction moyennée 


1 27 
Fiolr) = | F(r,6)d6. 
0 


T 
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sachant que 
27 ir ve : Ho re : 
cos!(8) sini+2(0)d0 = 0, si à impair où 7 impair, 
0 Fi si à pair et 7 pair, 


où à;; est une constante non nulle, on trouve que 
l i+j+1 
For) = 5 > Qij,20i5T : (3.15) 


où i et j sont paires. On voit que le polynôme Fig(r) a au plus [2] racines positives. 
De plus ce polynôme peut admettre exactement [2] racines positives simples, avec 


un bon choix des coefficients a;;. Cela complète la preuve du théorème (3.3.1). 


Exemple 3.3.1. Considérons le système 


à = y—e(V2 + 5x°?y + y°)y, 
Ts (3.16) 
ÿ=—t—e(l+ay+(8+ay — Sy )y). 
En coordonnées polaires (r,0) où x = rcosÛ, y = rsin0, r > 0, le système (3.16) 
devient 
F = —€ |V2r cos(0) sin(0) + sin(0) — r° cos*(@) sin(0) + r° cos(0) sin(0) 
1 
+ 4r* cos*(0) sin?(8) + a sin”(9) cos?(0) + r*sin?(0) cos(4) + r° sin?(4) cos(4) 
1 
+ 8r sin?(0) — - j sin?(8)| : 


ô=-1+ < — r sin(@) cos”*(9) — 3sin(0) cos(0) — 4r°* sin(0) cos* (0) — cr sin(@) cos*(4) 


+ Dr sin(8) cos(9) — SO 4 ar cos?(0) sin(8) — Vc0s2(8) + V3 + 12 cos4(4) 


— 2r? cos*(0) + r? , 


Pour déterminer les cycles limites, on résout l'équation 


3 
For) — 16” (—8 + 1] = (. (3.17) 
L'équation (3.17) possède une racine positive r — 2V2. D'après le théorème (3.3.1), 
le système (3.16) a exactement un unique cycle limite qui bifurque des orbites pério- 
diques du centre linéaire à = y, ÿ = —x. Voir la figure (3.1). 
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F 
" 
L 


CC 7 7 7 07 | 
mnt et dd | Ep 
ad a 8 gp Cp 


Th EE a a ve ve 4 4 à 
FE à à 4Ma X x % à 


“| , à 
R Ê + 

À | #& F Ÿ 
A 7 « CA 
XX ar Cr 
X % Fe / 
Le CA “ # 
Lo 7e AT # # 
K.* " éT Ca 
Ov vw VO PS 64 _ #— dE 


FIGURE 3.1 — Cycle limite du système (3.16) pour € = 0.0001. 


Preuve du Théorème (3.3.2). 


On utilise la méthode de moyennisation d’ordre deux, pour cela prenons f11, foi 
et ga les fonctions définies dans la partie précédente, de plus soit : 


fix, y) = D! Cast fa(x,y) = Ne. gax(x, y) sr. 
i+j=0 i+j=0 i+j=0 


Le système (3.12) en coordonnées polaires (r,0) où x = rcos6, y =rsin0, r > 0, 
devient 


— —E(A+eB), 
à 3.18 
0=—1--(4;+eB:), (3.18) 
a 
où 
> &jar +7 cos'(4)sin?*?(8 ns bjor ti cos'(0) sin/*1(0) 
VE i+3=0 
l 
+ D auarititt costt}(8)sini* (6), 
i+j=0 
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nm 
D Cij,27 TT cos! (0) sin/+?(9 #5 dir" cos’ (0) sin?*1(0) 


ne i+j=0 
+ 5 Cijar +51 cosi+1(9) sin/+1(0), 
i+j=0 


et 


n mm 
A — ÿ aj2r TT cos" *1(0) sin/*t1(0) + D. bjr T7 cos +t1(0) sin/(0) 
Le i+j=0 


- 5 aijir TT cos" (0) sin/*?(8), 


i+j=0 


nm mm 
Bi — D» C7 TT cos +1(0) sini+1(0) + ÿ dij2r T1 cos *!(0) sin’ (0) 
i+j=0 i+j=0 
l 
_ >» Cijar +1 cos’ (0) sin/*?(0). 
i+ÿ=0 
Considérons maintenant 0 comme la variable indépendante, on obtient 


dr 


Un ePi(r,0) + e*P(r, 0) + O(E), (3.19) 


avec 


Ce ()) — À, 
Pr, 0) = B — ! AA. 
T 


Déterminons la fonction moyennée suivante 


L, JF [la 
F. = — — F\(r, 0 ,0) + Pr, 0)| d6, 
= 7 [AGO + 0) 
où y1(r,0) = ie Fi(r,s)ds. Pour cela, on a besoin d’annuler la fonction moyennée 
Fo, donc on pose 4;;2 = 0 pour tout À pair et j pair. Maintenant, on commence par 
calculer 


d - ne , _ 
air, 0) = >» (à + 5 + Lajar 9 cos'(0) sin/T?(0) + ÿ (i+ jjbigart tit 
PEU i+j=0 


i odd or j odd 
l 


x cos’ (0) sin/T1(9) + ÿ: (+ j + 1aar +? cos **(9) sin/t”"(8), 
i+j=0 
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ensuite y1(r,0) qu’on écrit comme suit 


0 
nn) [| Finbat = y ++ 
0 
avec a 
yi(r,t) = | ÿ aij2r 71 Tcos"(t) sin/*?(t)dt 
0 5+5=0 
— &o27°(a110 sin(0) + @210 sin(30)) +. + ec DATES ne (@icie sin(0) 


+ Ace sin(30) + + ŒCaitat2H,, : sin((C1 + €1 + 2)0)) + ao1.27°(@01 
2 


€ 


T @201 cos(0) + @3o1 cos(30)) AUS ner OT (one + Aopiq cos(0) 
cos((p1 + qi + 2)0)) + anor (ani 
+ ou cos(20) + azu cos(40)) + +: + ar TETE + A2cim COS(20) 


T d3ciq cos(40) Fe RÉ cos((c1 QE T 2)0)), 


D 


où 

e c; est le plus grand impair et e, le plus grand pair tel que & + ei <n, 

e p, est le plus grand pair et q, est le plus grand impair tel que p1 + qi < n, 

° a; sont les constantes réelles obtenues durant le calcul de f} cos’(t) sin/*?(t)dt 
pour tout ? et 7. 


ge m 
yi(r,t) _ D bijor TJ cos'(t) sin/*1(t)dt 
0 


i+j=0 
= boo2 (&oo + &200 cos(#)) + bo2.,27? (ao + &9202 Cos(0) + 302 cos(34)) 
LP Lol Ce + pre2 cos(0) 4 Asp cos(30) 


cos((p2 + €2 + D0)) + bo1,27 (aid + 201 sin(20)) 


A(P2+e2 49) 


Pp2e2 


bo or Pr CT + Goprga (20) + Aspage SN (48) + + Àeratarts  g) 


X sin((p2 + q2 + D0)) + bio2T (au + &210 cos(26)) + b30.27° (ao + 230 cos(20) 
+330 cos(48)) Re bee2,2 Tire (aicses + Ü2c2e2 cos(20) DE U: çeateoti +1)c2e2 
X cos((c2 + €2 + D6)) + bi1,27? (a sin(0) D ne 911 sin(36)) + bi32r (aus sin(0) 


+ Ga sin(30) + às13 sin(50)) Lt bg 27 248 Ce sine) + oc Sin(30) 


+... + Gçeata2t2 qu sin((C2 + 2 + D0)). 
où 
e p2 est le plus grand pair et e, le plus grand pair tel que p2 + €2 < m, 
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e « est le plus grand impair et q est le plus grand impair tel que €2 + go < m, 
e «;;x sont les constantes réelles obtenues durant le calcul de 1 cos! (t) sin/*1(t)dt 
pour tout ? et 7. 


et 


CN 
y(r, t) ÿ aijar cos *t1(t) sin} (#)dt 
0 


i+j=0 
= â . Les+1 fx A 
= aoor (100 + 200 cos(26)) + + Gpges 17 PS TES (és + Àopzez COS(20) 


A 


Œrstest? 
2 


+1)p3e3 


p3+q3+1{ A i À i À i À 
dsail Ce sin(0) + Gopags Sin(30) + Àspags SIN(50) + … + À(patasts) 


cos((p3 + e3 + 2)0)) + ao117? (ésor sin(0) + &2o1 sin (30)) 


se p3q3 
X sin((p3 + g3 + 2)0)) + &10,17? (äuo + 210 cos(0) + 310 cos(34)) 


ca+e3+1{ 


+... + Gcses 17 (ércses + Àcses CO8(0) + Àscses COS(30) + … + Àçestes+ 


X cos((c3 + €3 + 2)0)) + &1117° (and + Ào11 sin(40)) + &1317° Cr + À 213 sin(20) 


+2)c3e3 


+313 sin(40) + À 413 sin(60)) SET PT LE | sea 0 + ous sin(20) 


+ Ü3cags SIN (40) + + À(cy+as)esgs SIN((C3 + 3 + 2)0)). 

où 

e p3 est le plus grand pair et e3 le plus grand pair tel que p3 + e3 < l, 

e c; est le plus grand impair et q3 est le plus grand impair tel que €3 + q3 < L, 

e «;;, sont les constantes réelles obtenues durant le calcul de Le cosi+1(+t) sin} (#)dt 
pour tout 2? et 7. 


Finalement 
yi(r, 0) = a10 27? (ao sin(0) + @10 sin(30)) + + Que 27 TA Ce sin() 


+Ocies SN(30)+... + Qatar, Sin((c + e1 + 2)0)) + @01,27? (ao + @201 COS(0) 
2 


+@301 cos(34)) + + Qpig 27 PTT Cr +O2piq COS(O) + Aapiu COS(30) 


Qta+2+3 cos((p1 + qi + 2)6)) +a1127 Cr + @211 COS(20) + a 311 cos(48)) 


— p1q1 


+q+1 
ag ar tat (arc 2e COS(20) + A3c1m COS(48) + … + data +20 


X cos((c1 + gi + 2)6)) + boo,2 (a00 + 200 cos(#)) + bo2.27? CT + 202 cos(0) 
+ &302 cos(34)) +... + be,.2rP2762 (ose + Gpye COS(0) + Aapye> COS(30) 


De) 


cos((p2 + €2 + D6)) + bo1,27 (aid + 201 sin(20)) 


p2e2 


bag o7P2+8 CT + Gapage Sn (20) + Gpaga SIN(40) + + Açraterts 
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X sin((p2 + q2 + D0)) + bio2T (au + 210 cos(26)) + b30.27° (ao + 230 cos(20) 


+1)c2e2 
X cos((ca + €2 + D6)) + bi1,27° (a sin(0) Sr 911 sin(30)) + bi3,27* (aus sin(0) 
+äns sin (30) + ân3 sin(50)) ++ bag ar 248 (és sin(8) + À2cvç Sin (30) 


+330 cos(48)) us bises re2+e2 (aicses + Ücpe) cos(26) haie eteoti 


FF ete), £ sin((c2+42+1)0)) + @00,17 (d100-+ 200 cos(26)) ++ @pges 1723834 


2q 


x (éme + Âopses COS(20) + … +Q(ratest2 cos((p3+e3 +2)0)) + ao1.17? (ér sin(0) 


+1)psez 
+201 sin(30)) + + Qpsgs 1773781 Ce sin(0) + Gopsgs SN (30) + Àspsgs SM(50) 


A 


Arstasts) eq sin((p3 + g3 + 2)6)) + &10,17? (äro + 210 cos(0) + 310 cos(34)) 


c3+te3+l1l{ A A À 
Ages ar 343 (ses + Àocses CO8(0) + Àscses COS(30) + … + Âçeatesti 49) eg 


X cos((c3 + €3 + 2)0)) + &1117° (and + 11 sin(40)) + &1317° (ar0 + À 213 sin(20) 
+313 sin(40) + À 413 sin(60)) an LE | disst + Ses sin(20) 


Agen SN (40) + + Éfeuranjeng SIN (C3 + @3 + 2)0)). 


En utilisant les résultats des intégrales de l’Appendice, et notons 


Hi(r) = - u LR (r,Ojy(r, o)] db, (3.20) 


on trouve 


Hi(r) = :| Duras (G+7+lagert laio2r? (aoAt ca anoAÿ,) 


à pair et j impair 
citei+l 1 3 c1+e1+2 
TT os + Qcie:,27 L L (arcses A + Acier A}; + des + ACatatatte,e, Ai ) 
mn ; : i+j 2 1 3 
En Dies (i nn De 1a;or ' Laoi.2r (18; ou 1) 
i impair et j pair 


+qi+l 1 3 p1+91+2 
So Du d'A (ame 83 + Ascpiu Pi + + + Aatarass, Bi 


n : : i+ & D1 2{ D1 = 3 
LE D GT Da sr [bo02 7) + Do2,27 Cr LE 7) 
i impair et j pair 
ve p2+e2 [ à D1 x 3 = p2+e2+1l 

LT bar (ame: B} + @302B;; es Gçrate2 + 2)ppen Pis 
n ; 5 it = 2 
+ re (sbiasssr te Lbor.2r (rois NC GriCà) 
i impair et j impair 
p2+g{; e = 2 « 4 À ; p2+q2+1 

ce + Opoge.2T (éarae + GopagOij + Aspago Ci + ee + rater), 0, Ci 

n es +5 2{( x... A1 Fe A3 
+ De (+3 +1la;ort Lonar (anÂt + âAÿ,) 


i pair et j impair 
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c2q2 1] 


Dre (+ j + Laiart Lao ar? Cr + Gao1 À? ) 


i pair et j impair 


p3+q3+1{( A Â1 2 13 à Àp3+aq3+2 
4 Apsqs 17 (Give À5 ua Àpags Àÿ; Sr n Gcéstasts) A 


ES A1 ES A c2+g2+1 
«+ PTT ARE CerE + G2c2g À; Ra çcotot2), A°° )| 


paq3° 1j 
D (+7 + ljajart laroar? (ab; + 310 B2 :) 
i impair et j pair 


c3+e3+1 D C3 +e3 +2 
nn LA (re B} + Éscses PB? +... + Âçestesti  ojes Pis 


nm | 3 4 
Danse G+i+lja;er do fauar C7 S: émC) 
i impair et j impair 


r3+a3+1 AC3+93+2 
+ + Gcsgs,17 7. Ce F Gaesgs C2 AD ur À(c3+g3) csas Ci 


m LL i+j—1 
Br (+ j)bijar 9 laro2r? (auoD + @210DÀ ;) 
i pair et j pair 
ci +e1+1 cit+ei+2 
DS dre ol 1+€1 (ace D} + Ocies D; Hat Giécererer, _. D; ) 
m Là i+j-1 2 
er (à LE J)bijor d Laoi2r (aE} SF as E$) 
i impair et j impair 
+qi+1 1 3 pi +g1+2 
+ Gprgi 2 PTT (ame Ei 4 pig E5; Ar Awtatts, E;; 
m Fe +j—1 3 
re (à = JJbijortt3 lauar (on air ou F$ “5 Q311F% +) 
i pair et j impair 
ci +q1+1 4 C1+g1+2 
+ Gciq1,27 1+q1 CT + Ganrs + A3c1q F;; see ice ii ne 
m Hi es i+j—1 = pl 3 
Do (+ j)bgort 1 bon ( &200É4, ) + +bo22r? | &202É 4 + &302E 5, 
+3 J 4 J 
i impair et j impair 
p2+e2 [p2+e2+1 
+ Dpses, 27 Cr + Gares EE +. + à» 22722 rer 


sn. (à je Dbsrt a Lborar (ao; + 1) 


i impair et j pair 


a = p2+q2+1\ | 
Het bpaga,2r 72 te Cr + Gapaqo C3; TR etats), (LE ) 


ve (è + jbisart tit (bo, 27 (ao: + à10 F2 à 


i pair et j impair 


c2+e2 prc2+e2+1 
She be 27° Cr? + än0F2 DE OR erteoti ours 


m : : i+j—1 21. 1 
Es (è + j)bijar 9 Loar (and, + Ga D},) 
i pair et j pair 


c2+q{f DL = nc2+a2+1l 
. + bar C2 + De ET «+ çatete on Di 


m ; S i+j—1 2 
Dire + )ber Lao, ir (008: + G00 F; 
i pair et j impair 


ps+ez+1/f d 2 pp3+e3+2 
: + Apses 17 (&ipses A ci Éopses FE ++ Àçrstest2 jme Là 


P2q2 


60 


3.3. BIFURCATION DES CYCLES LIMITES DES ORBITES PÉRIODIQUES 
D'UN CENTRE LINÉAIRE PERTURBÉ PAR UNE CLASSE GÉNÉRALISÉE 
D'ÉQUATION DIFFÉRENTIELLE DE LIÉNARD 


mr HR. i+j—1 2/4... Di s.  D3 
D (é + j}bior ts faoiar C2 + Ga D) 


i pair et j pair 


p3+q3+1{ 4 1 À 3 2 p3+93+2 
+ À Gpsgs,17 De + Àopsgs Di + + Àçratests Di 
m een i+j—1 2/4. 1 s . E3 
Durs (C+J)bar? faror (aoË, + éoËt,) 
i impair et j impair 
c3te3+1f A pl 2 3 A cs +e3+2 
+ Acses,17 Âacses E = Àscses Ei Re Far Hosts 


m - : i+j—1 3{/: 2 A4 
Die (i ua J)bijar “ fanrar (ao Le éG) 
i impair et j pair 


c3+g3+1/{ z A 72 2 4 2 ca+g3+2 
ce À deg 17 CPC £E Â2cags Gi F Âcsgs Gi LE À(c3+98)c398 


ij 
= CE DTA LE (nor? (onoAy LE anoA) 


i impair et j pair 
citei+l 1 3 ci+e1+2 
… + Qcie:,27 1 L (are: A} + ce H}; + …. + Qate+4,. Hi 
2 


l : : i+3 2 1 3 
Dur (+7 +1lja;irt Laoi.2r Co1l;; + Q3o11;; 
i pair et j impair 


+q1+1 1 3 p1+q1+2 
LR Gp PET CE SA pra Li ++ Qeitn+2+3,, on li 
2) 


l : : i+j 3 2 4 
Dr (i+ 5 + Las Jaurar Cr + Qi K}; + a311À;; 
i impair et j impair 


c+qitl K€ +q1+2 | 


c1q1 ij 


De (+5 +Lairti [b002 (ä200/3) + b2,27° (ay LE 15) 


i pair et j impair 


2 
an ei q,27 CRC + ion À; + + Œ(cy+a1+2)+2 
2 


pates [| % 1 = Tp2+e2+1 
Sr Dp2e2,27 CE ++ dote +2)pae Lij 


nus (à +3 + L'air Lborar (érour, + Go Là,) 


i pair et j pair 


[L?2 +g2+1 
p2q2 1] 


ee (+3 + Dagart lbroar (aoû + Gao?) 


i impair et j impair 


a = 2 + sr c2+e2+1 
+ baser en Cr + Ace K}; ++ Gçeatertt | 1jone Ni ) 


l ne Se +5 2/$.. 11%. D 
Den (ri bDaar te Lonar (an ñy + Go A}, 
i impair et j pair 


LR . Cri z crc2+q2+1l 
SA +- baser 2792 CRE + .. + çotete) oo Hi 


p2+q{A sf 2 À 
ce À Opag2.27 Cr + Gopago Li + + CHE" 


l : : à . k Re A 2 
Des CT Last Laooar (aoû: Æ Go 8 ) 
i impair et j impair 
pates+1f SAS F2 = -p3-+e3 +2 
+ + Gpses il C7 À. pes Ki Has Àçratest2 mes Pi 


l si i+j 2fx T1 a T3 
ne, GET Laorar @o1 À; + Goo A}; 


i impair et j pair 
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i pair et j impair 


+l)cses 19 


i pair et j pair 


LL c3+g3+1 
au Ac3g3, 17 CET or Gcsgs L2 Sun À(c3+g3) 


Maintenant, on calcule 


nm 
Ü}= > C7 TT cos’ (0) 0) sin/*?(9 3% d; 


i+j=0 i+j=0 


ps+g3+1{ CT 1 À T3 à 
+ À Apsqs 17 (éme À 5 Âpsgs Hi; Le A(rstasts 


; p3+q3+2 
43 pe D )| 


l pe +5 2/5... F1 SU 
JL Dee (à +9 + Dossier Laroar (ao + à3101;; 
À T 2 
+ + Gegeg ar 8781 (ref + - À çestestt Le tes+ )| 


NSP GED lauars Cr + énlé) 


T c3+q3+2 
cags Li ) | L 


27 cos! (0) sin/*1(8) 


dE S Ci, rl os°+1(0) sin*1(9) 
i+j=0 
n n 
EH > ) Gad TRE cost T1 (0) sin *+3(0) 
i+j=0 k+h=0 
i impair ou j impair k impair ou h impair 
n m 
_ ) > Qij20kn,27 ITA cos TF+1(0) sini+#+2(9) 
i+j=0 k+h=0 
i impair ou j impair 
n l 
— Ù > Qij2an ir TI TETREL cositE (0) sini+#+4(9) 
i+j=0 k+h=0 
i impair ou j impair 
m n 
_ > > bijaagn.or TITEER cositEt1(0) sini+h+2() 
i+j—=0 k+h=0 
k impair ou À impair 
| ; bobinor FTREl cos T#+1(9) sin t#+1 (8) 
i+3=0 k+h=0 
+ Ù Ù bijaarnar TER cos TE (0) sini+*#+3(9) 
i+3=0 k+h=0 
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l n 


_ ) > aijaagh or FEAT cost #t2(9) sin +2(0) 


i+j—=0 k+h=0 
k impair ou h impair 


l m 
— ù ) aijabrn or TER cost +2(0) sini+#+1(0) 


i+3=0 k+h=0 
l l 
+ +k+h+I i+k+1 DE] 
da Ù Ù AN NS ne) ut 
i+3=0 k+h=0 
et notons 
1 27 
Br) = — FP2(r,0)d0, 
27 Jo 
on obtient 
1 n m 
Es i+j+1 it 
Hi(r) = 9 | ) Cia Qi + ) dij27 Bij 
i+j=0 i+j=1 
à pair et j pair i pair et j impair 
l 
+ ) Car 6, 
i+j=2 
i impair et j impair 
n n 
+j+k+h+1 
_ ) ) Qij2@kh,27 7 Mijkh 
i+j=1 k+h=1 
i impair et j pair k pair et À impair 
n n 
i+5+k+h+1 7 
— > > Qij2@kh,27 7 Men 
i+j=1 k+h=1 
à pair et j impair k impair et hk pair 
n m 
+5+k+h 
= > > Qij,20kh,27" 77 Nijrn 
i+j=1 k+h=0 
i impair et j pair k pair et À pair 
n m 
i+j+k+h N\ 
Du ) ) Gij,20kn,27 Ë Nijuh 
i+j=2 k+h=1 


i impair et j impair k pair et hk impair 
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3.3. BIFURCATION DES C 
D'UN CENTRE LINÉAIRE 


YCLES LIMITES DES ORBITES PÉRIODIQUES 
PERTURBÉ PAR UNE CLASSE GÉNÉRALISÉE 
D'ÉQUATION DIFFÉRENTIELLE DE LIÉNARD 


n m 
i+j+k+h À 
CY ) ) Qij20kh,27 ï Nijuh 
i+j=1 k+h=—2 
à pair et j impair k impair et hk impair 
n l 
à+j+k+h+1 
+ ) ) Géo AT Fete 
i+j=1 k+h=1 
i impair et j pair k impair et h pair 
n l 
i+j+k+h+1 D 
+ > > Qij2Qkh17 7 FRE 
i+j=1 k+h=1 
à pair et j impair k pair et À impair 
n l 
i+j+k+h+1 Rp 
+ > > Gij2@kh 177 Fotr 
i+j=2 k+h=—2 
i impair et j impair k impair et hk impair 
m n 
i+5+k+h N 
Te ) ) bij2@xn,27 : Nijuh 
i+j=2 k+h=1 
i impair et j impair k pair et hk impair 
m n 
i+j+k+h 
LE ) ) bij24gh,27 s Ni 
i+j=0 k+h=1 
à pair et j pair k impair et h pair 
m n 
i+j+k+h À 
= ) ) bij2@xn,27 d Nijuh 
i+j=1 k+h=—2 
à pair et j impair k impair et hk impair 
m m 
i+j+k+h-1 
LL ) ) bij 20kn.27 : Qibh 
i+j=0 k+h=2 
à pair et j pair k impair et h impair 
m m 
i+j+k+h-1 
ca ) ) bij 20kn.27 ÿ ur 
i+j=2 k+h=0 


i impair et j impair 


m 


> 


i+j=1 
i impair et j pair 


k pair et h pair 


m 


> 


k+h=1 
k pair et À impair 


i+j+k+h-1 
bij,20kn,27 Que 
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3.3. BIFURCATION DES CYCLES LIMITES DES ORBITES PÉRIODIQUES 
D'UN CENTRE LINÉAIRE PERTURBÉ PAR UNE CLASSE GÉNÉRALISÉE 
D'ÉQUATION DIFFÉRENTIELLE DE LIÉNARD 


m m 
i+j+k+h-1/ 
LE ) ) bi;20gn,27 . Oh 
i+j=1 k+h=1 
à pair et j impair k impair et h pair 
m l 
i+j+k+h 
+ ) ) bij2agna17 7 Rth 
i+j=0 k+h=1 
à pair et j pair k pair et hk impair 
m l 
i+j+k+h 
+ ) > bisopnart Rijkh 
i+j=1 k+h=0 
à pair et j impair k pair et h pair 
m l 
i+j+k+h p 
+ ) > bij2axn17 7 RER 
i+j=1 k+h=—2 


iimpair et j pair k impair et h impair 


l m 
i+5+k+hT) 
_ ) ) Qij10kh,27 : Dh 
i+j=1 k+h=—2 
iimpair et j pair k impair et À impair 
l m 
i+5+k+hT) 
| ) ù Qij10kh,27" 77 Uijkh 
i+j=2 k+h=1 
i impair et j impair k impair et À pair 
l l 
à+j+k+h+1 
+ ) > LOT Vijkh 
i+j=0 k+h=2 
à pair et j pair k impair et hk impair 
l l 
i+j+k+h+1T 
+ ) ) RU 170 LINE Vijkh 
i+j=1 k+h=1 
à pair et j impair k impair et h pair 
l l 
i+j+k+h+1T 
Fe ) ) Qij1akhA17 7 Vijkn |. 
i+j=1 k+h=1 


iimpair et j pair k pair et hk impair 


Finalement, on remarque que H\(r) + H(r) est un polynôme en r? de la forme 
suivante. 


Hair) + ar) = Pi(r?) + Pr) + rt Por) | 
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3.3. BIFURCATION DES CYCLES LIMITES DES ORBITES PÉRIODIQUES 
D'UN CENTRE LINÉAIRE PERTURBÉ PAR UNE CLASSE GÉNÉRALISÉE 
D'ÉQUATION DIFFÉRENTIELLE DE LIÉNARD 


P,(r?) est un polynôme en la variable r? de degré 
O(m) + E(t) — 1, EC) +O(i) -1 
2 l 2 


PAU +0) — LE 


P;(r?) est un polynôme en la variable r? de degré 


À = max{| Re AL 


+ On) —3, Om) + Efn) —3, On) +O(L) —2 
PUR) $, [OUR + EU LE, OU +S 


Fm) — 2, EC) ee = 3, [O(m) 0 3 


P3(r?) est un polynôme en la variable r? de degré 


À2 = max{O(n)-1, EAU l,{—1, 


E(n) + E() — 4 
2 ]; 


ef 


avec 

e O(i) le plus grand impair inférieur ou égal a 1. 
e Eli) le plus grand pair inférieur ou égal a 1. 

e |.] désigne la partie entière. 


Donc 


1 
Fo(r) = 9 


ns 
T 


LP?) + r2P2(r2) + ré Pa(r?)]. 


Pour trouver les racines positives réelles de F39, nous devons résoudre une équation 
algébrique en r? de degré À = max{ A, À + 1, À3 +2}. Cela donne que F2 a au plus 
À racines positives. Cela complète la preuve. 


De plus, nous pouvons choisir les Coefficients à; 1, @;,2, bia, Ci 1, C2, di; 2 de telle 
façon que F9 admet exactement À racines positives. 
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3.3. BIFURCATION DES CYCLES LIMITES DES ORBITES PÉRIODIQUES 
D'UN CENTRE LINÉAIRE PERTURBÉ PAR UNE CLASSE GÉNÉRALISÉE 
D'ÉQUATION DIFFÉRENTIELLE DE LIÉNARD 


Exemple 3.3.2. On considère le système 


ER R 2x? — y — E(—2xy — y°)y, (3.21) 


ÿ = —2 — e(4ry + (—2* — y) — (a — y + (ay + y? — z)y), 


En coordonnées polaires (r,0) où x = rcosô, y = rsin0, r > 0, le système 


(3.21) devient 


— 1 — cos(8)sin(6)r + (sin(@) — 4sin2(9) cos(8)}r? + 2 sin(0) cos” (8)r | 
+ &2[sin2(8}r + (—sin(8) cos?(8) + sin?(8) cos(8)}r* 

+ (2c0s(8) sin(8) + 3c02(9) sin?(9) — 2 cos? (4) sin(4) — sin?(0))r*} 
is. É — cos2(0) — 4sin(8) cos? r + (cos() sin(8) (3.22) 
+ 2cos4(8) — 2 cos?(8))r : 4e? | cos(0) sin(0) + (— cos#(6) 

+ sin(9) cos?(8))r + (3 sin(0) cos?(9) — 3 cos(0) sin(0) — 2 cos(8) 
+ 3 cos?(8) — 1yr?| 


Considérons maintenant 0 comme la variable indépendante, on obtient 


_- = &F\(r,0) + Pr, 0) + O(e5), (3.23) 


où 


Fi(r,0) = —2r* sin(@) cos*(@) + cos(6)r sin(9) — r? sin(0) + 4r? sin?(0) cos(0), 
P(r,0) = —6r“sin(@) cos“ (0) + 20r° sin(4) cos”(9) — 2 sin(@)r° cos(@) + 6r“sin(@) cos*(8) 
— cos*(@)r sin(0) — 4r° sin(0) cos/(9) + 4r° sin(@) cos” (0) — 18r*sin(0) cos°(0) 
+ 2sin(0)r° cos?(0) + cos()r sin(0) — r?sin(0) — 8r* sin?(0) cos*(4) 
+ 3r? sin?(0) cos(0) — r sin?(9) — 8 cos*(0)r” sin?(0) + r° sin?(0) + 2r* sin?(0) cos” (0) 
+ 16r“ sin°(@) cos”(0) — r sin?(@) cos(0) — 2r° sin?(8) cos*(0). 
1 


Notons d’abord que : Fio(r) = 5x Fee 
T 
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3.3. BIFURCATION DES CYCLES LIMITES DES ORBITES PÉRIODIQUES 
D'UN CENTRE LINÉAIRE PERTURBÉ PAR UNE CLASSE GÉNÉRALISÉE 
D'ÉQUATION DIFFÉRENTIELLE DE LIÉNARD 


Pour déterminer les cycles limites, on résout l'équation 


nue ef Do. [Ress + F0 do 


2m à . (3.24) 
= RAT — ri — 5) —=Ù 


1 1 
L'équation (3.24) possède deux racines positives r1 = 5 V10 + V2 
1 1 
et Ta = 5 V10 _ 5 V2. D'après le théorème (3.3.2), le système (3.21) a exacte- 


ment deux cycles limites qui bifurquent des orbites périodiques du centre linéaire 


à = y, y = —x, en utilisant la méthode de moyennisation d'ordre deux. Voir la 
figure (3.2). 
A AAA AN A 7 “a a+ 
4 AMAAAN AN "D'R" 
# 1 1 : : L | \ 
; Ÿ A #1 4 _ x V \ 
# AA A 5 7 VO \ 
daté « X x 
4 1 A1 À 1 8 9 LA 
l 4 s 4 à À 
A1 14 9% VU A 
ae D \ X à 
1 411 1 1 Ps \ À \ 
PE. OR 
4 4 4 4 74 NT 
RP L Ÿ YX Ÿ 
4 4 fs 1 | 
DES) RE Y + + 
4 4 4 À | | | 


La: 


“| AT LEE 
à Po EDR 
a FT PRET D Te 
 * ERFÉFÉFF 
kKhA TE ÉÉÉÉY 
X* CNE HE 
LR DEN NES 
CE CCS ll 
> % CN ddr 
x % Ca LE 8 a pp 


FIGURE 3.2 — Les deux cycles limites du système (3.21) pour £ = 0.0001. 
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- Chapitre 4 : . 
Bifurcation des cycles limites d’un 


anneau périodique pour un 
système différentiel plan 


Dans ce chapitre, nous allons appliquer la méthode de moyennisation proposée en 
2004, par Buicà et Llibre [6] pour étudier le nombre de cycles limites qui bifurquent 
des trajectoires périodiques d’un anneau périodique d’un système différentiel plan. 
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4.1. QUELQUES DÉFINITIONS 


4.1 Quelques définitions 
On considère le système différentiel 
(4.1) 


Soit ( un ouvert du plan de phase et (x(t), y(t)) une solution (ou une trajectoire). 
Notons par A l’ensemble 


A = {te R/(x(t),y(t)) € Q. 


Définition 4.1.1. (Intégrale première). L'application H : Q — R est appelée in- 
tégrale première du système (4.1) sur Q si elle est constante sur les courbes solutions 
(æ(t),y(t)) du système (4.1) contenues dans Q c-à-d 


H{x(t),y(t)) = cste,Vt € nr 


Définition 4.1.2. Une fonction différentiable (x,y) — H(x,y) est une intégrale 
première du système différentiel (4.1) si : 


OH OH 


Remarque 4.1.1. On dit que le système (4.1) est intégrable sur un ouvert Q s'il 
admet une intégrale première sur (1. 


Exemple 4.1.1. Soit le système différentiel 
à = ax — bxy 
ÿ = cry — dy 
définie sur U =]0; +co[x]0; +oo[. La fonction 


F:(x,y) EU = cx + by - dinx - alny 


est une intégrale première du système. En effet 


D'où le résultat. 


70 


42. THÉORÈMES 


Définition 4.1.3. (facteur intégrant). On appelle facteur intégrant du système 
(4.1) sur Q associé à une intégrale première H une fonction de classe C* 
R:Q—R non identiquement nulle telle que 


O(Rf) | O(Rg) 
Ox ï Oy 


=D édo(RX)= 0. 


Définition 4.1.4. (Centre isochrone). Un point singulier est un centre si, dans 
un certain voisinage de ce point toutes les orbites sont fermées. 

Un centre est isochrone si le temps de parcours de ces orbites fermées est toujours le 
même. 


4.2 Théorèmes 


On considère le système différentiel autonôme plan 


ne 


ÿ = Q(x,y), (42) 


où P,Q : R? — R des fonctions continues. Supposons la condition (A1) suivante : 
(A1) : le système (1.2) a un anneau périodique autour du point d'équilibre (0, 0), 


T,:{(x,y) ER: H(x,y)=h, hk<h<h,}. 


où H est une intégrale première du système. 
Désignons par u = u(x,y) le facteur intégrant pour le système (4.2) correspondant 
à l'intégrale première AH. 


On perturbe le système (4.2) comme suit 


{ à = P(x,y) +ep(x,ye), (4.3) 


y = Q(x, y) T Eq(x, y, 2), 


où p,q:RXR—R 
En 2004, Buicà et Llibre [6] ont proposé une manière pour appliquer la méthode de 
moyennisation dont le but est d’étudier le nombre de cycles limites qui bifurquent 
des trajectoires périodiques du système pour € = 0. 


La première étape est d'écrire le système (4.3) sous la forme standard pour ap- 
pliquer la méthode de moyennisation, i.e. sous la forme 


& =ef(t,r) + q(f, x, €). 
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42. THÉORÈMES 


Le système différentiel (4.3) écrit avec cette forme standard décrit la dépendance 
entre la racines carrée de l'énergie À = V/h et l’angle w des coordonnées polaires. 


Théorème 4.2.1. Supposons que (Al) est vérifiée pour le système (4.2) et que 
rQ(x, y) — yP(x,y) À 0 (4.4) 


pour tout (x,y) dans l’anneau périodique. 
Soit p : (Vhe, Vh:) X [0,27) — [0, oo) une fonction continue telle que 


H(p(R, &) cosw, p(R,v)sing) = F?, (4.5) 


pour tout R € (he, Vh:) et tout w € [0,2r). 


Alors, l’équation différentielle qui décrit la dépendance entre la racine carrée de 


l'énergie R = Vh et l’angle $ pour le système (1.3) est 


dR _ u(æ +y)(Qp — Po) 
d 2R(Qx — Py) + 2Re(qx — py) 


(4.6) 
où x = p(R,v)cosvy et y = p(R,v)sine. 


Preuve 


D'’aprés la définition de l’intégrale première et du facteur intégrant, nous avons : 
OH OH OH 0H 
por = P = 


y — 1, 5x —HQ. 


On définit la fonction 

G(r, R,4) = H(r cosy,rsing) — R?, 
en tout point (r,) de l'anneau périodique et pour tout R € (he, Vh:). 
On a 


0G OH 
Ôr x 


COS D + 


Où ZT = Tr Cosp et y = rTSIN. 
Pour tout (ro, 0) dans l’anneau périodique, il existe un À tel que G(ro, Ro, vo) = 0. 
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42. THÉORÈMES 


0G 
L'hypothèse (1.4 =. (ro, À 0. 
ypothèse assure que De (ro; Ro, Po) 7 


D’aprés le théorème des fonctions implicites, au voisinage de chaque point (Ro, @o) 
il existe une unique fonction continue p = p(R,) telle que l'équation (4.5) soit vé- 
rifiée. 


Donc cette fonction est bien définie sur (he, Vh.) x [0,2T) et satisfait (4.5). 


La relation entre la racine carrée de l'énergie et le temps est donnée par 


R(E) = V H(x(t), y(t)) 


et entre l’angle & et le temps est 


pour tout (x(t),y(t)) e lr,tEeR. 


On obtient alors 


p--M@r— Pa), _(Qx—yP)+e{ax — py) 
TRE PCR ; 


éliminant le temps dans les deux expressions, on obtient l’équation (4.6). 
P P 


La démonstration du théorème qui suit est une conséquence du résultat du théo- 
rème 4.2.1. 


Théorème 4.2.2. la fonction f : (Vhe, Vhs) x (—-€/,€r) — R définie par 
jte es F(t,x(t,z,e),e)dt pour l'équation (1.6) est donnée par l’expression sui- 
vante 


__ ff u(x?+y)(Qp- Pa) 
J(RE) = cf 2R(Qx — Py) + 2Re(qx — py) " 


et la fonction fi : (Vhex, Vhsx) — R définie par f1(z) = Un F\(s,z)ds pour l'équation 


(4.5) est 


(4.7) 


A (R) = [ (x? +y?)(Qp— Pa), 


2R(Qr-Py (48) 


où pu = u(x,y) est un facteur intégrant pour le système (4.2) correspondant à l’inté- 
grale première H, avec x = p(R,v)cosv et y = p(R,v)sine. 
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4.3. BIFURCATION DES CYCLES LIMITES D'UN CENTRE ISOCHRONE 
VIA LA MÉTHODE DE MOYENNISATION 


4.3  Bifurcation des cycles limites d’un centre iso- 
chrone via la méthode de moyennisation 


Comme il a était présenté dans [6], soit le système 


. 

T—=—-Y+T, 

_ (4.9) 
Y = X + Ty, 

qui possède un centre isochrone à l’origine et une intégrale première 


x? + 
H(s;u) = ;- Notons que 4 — 0, h; = 1, et la fonction p satisfaisant les 


(1+y) 


R 
hypothèses du théorème 4.2.1 est p(R,@) — TE Rsno Pour tout 0 < R < 1et 
— Rsing 


p € [0,2x). 


On perturbe le système (4.9) comme suit 
ue 
ÿ 


pr, y) = ar — a3x° + (2a2 + a5)ty + a6y° 


—V no & no ep(x, y), 
x + œy +eEq(e, y), 


(4.10) 


où 


et 
q(x, y) = ay + Q2t° + Q4Ty — aY°. 


L’équation (4.6) correspondante à ce système est 


dR = Qui + a(p)R? + b(o)R° (11) 
dé  1—Rsing+ec(v)R | 


où 


a(p) = (—2a1 + 3a2 + as) sin @ + (aa + ag) cos g + (—5a2 — as) sin° w 
+ (—as — d4 — 06) cos o, 
bg) = a1 + a2 + (—a1 — 2a2) cos” © — ai cos psin y, 


c(y) = (az + a4) sin ç + (—3a2 — as) cos ç + (—a3 — a4 — a) sin° 
+ (4a2 + a5) COS” . 
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4.3. BIFURCATION DES CYCLES LIMITES D'UN CENTRE ISOCHRONE 
VIA LA MÉTHODE DE MOYENNISATION 


L’équation (4.11) s'écrit sous la forme 


dR 
a — eF(, R) + = G(p, h,ë); 
dp 
avec R )R2 b( }R# 
ahR+a + 
F(R, p) = : é Le , 


1—-Rsny 


lai R + a(p)R? +b(6)RIc(e)R 
(1 — Rsinv)(1 — Rsiny+ec(v)R) 


G(R, P; €) on; 


c’est à dire sous la forme standard de la méthode de moyennisation d'ordre un. 


Pour appliquer le Théorème de moyennisation d’ordre un, on a besoïin de la fonction 
hit) = 4 Fi(s,2)ds, qui est pour notre problème fi :[0,1[- R, 


nl 


2 Q12 + a(o)2? + b(w)z 
, — zsinçg 


On calcule cette intégrale par maple, on obtient 


1 
he EE LS + (Gaz + as — 2a1)2° V1 — 22 — (10a2 + 2a5)2? 


— (2a5 + 8a2)V 1 — 22 + 8a2 + 245]. 


En considérant la nouvelle variabla £ €]0, 1[ définie par z = 4/1 — £?, on obtient 


fi(v 1 — €?) — xt ci £)(2asé? + (2a1 — 4aa _ as)Ë + 24; + 2@) + as). 


Notons que z € (0,1) est une racine de f, si et seulement si € €]0, 1[ est une racine 
du polynôme g(£) = 2a2£? + c1Ë + co, où 1 = 2a1 — da — a5 et Co = 2a1 + 209 + 45. 
La fonction g possède au plus deux racines dans intervalle ]0, 1{, ceci signifie que la 
fonction f. admet au plus deux racines. 
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4.3. BIFURCATION DES CYCLES LIMITES D'UN CENTRE ISOCHRONE 
VIA LA MÉTHODE DE MOYENNISATION 


PAR æem EE = — — T = - 
/ ” Ov % + 
EH # 1 
/ % + à 
EE # L | 
f y % D % À» 
ÉÉ EF 
{ f x ® NX % 
FÉF | 
1 | / ù Où D à 
ÉÉF F \ \ à 
1 | / R R Rh R 
F F É À SR 
{ | / h Oh & à 
FFF PF VF à 
| / 5 D à à 
CRC AE \ 
| | | 4 4 4 à 
4 l 1 1 l 
Il | g': 
#4 % +3 | 
\ | 4 4 4 4 
Ÿ Ÿ % (| } (| 
\ 4 À 4 4 4 4 
SOS, | ' SR OT ART 
\ À 4 4 4 4 
Ÿ % À ÉTÉ 1 
\ À # A À A 
\ à à F ff) 
VU f #1 1 4 
X 4 *x 4 / / 
\ À f 4 4 4 
\ à « 
\ 1 1 A A 
k à x 
4 + 
AS" 1 À A 


FIGURE 4.1 - Deux cycles limites du système (4.9). 


Donc au plus deux cycles limites peuvent bifurquent de l’anneau périodique du sys- 


tème (29). 
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Conclusion et perspectives 


Cette thèse de doctorat porte sur un aspect important de la théorie qualitative 
des systèmes différentiels planaires, à savoir les cycles limites. 


L'importance de déterminer le nombre de cycles limites fait l’objet de la deuxième 
partie du 16°" problème de Hilbert. 


L'application de la méthode de moyennisation pour étudier les cycles limites des 
systèmes différentiels dépendant d’un paramètre, s’avère ces dernières années trés 
fructueuse, elle à donné de bon resultats et à été appliqué par différents Mathéma- 


ticiens. Notons que nos différents calculs ont été réalisé par le logiciel Maple. 


Nous continuons à travailler sur des problèmes analogues. 


ré 


Appendice 


Dans cet appendice, nous définissons certaines formules écrites dans la preuve du 
théorème 3.3.1 et dans celle du théorème 3.3.2, pour à > 0 et j > 0. 


27 _ | 
| cos’ (0) sin/*?(0)d0 — { He RADAR EL DES 
l 0 si non, 


27 
| cos’ (0) sin/*?(0) cos((2h)0)d0 = 0 sii impair ou j impair, À = 0,1,2, 
0 


| cos (9) sin (8)48 . FPS SR 
0 
0 


si non, 


sii impair et 7 impair, 
si non, 


27 
| cos! (9) sin/+2(0) x 0d0 = { “ 
0 


sii pair et } pair, 
si non, 


cos'*1(0)sin/*t1(0) x O0d0 = { se 


AOij Si? impair et j pair, 
0 si non, 


cos’ (0) sin*1(9) x 0d0 = { 


27 sb 5 RES : 
| cos/*1(0) sin/*1(0)d0 — { Rôij Sii impair et j impair, 
0 0 si non, 


te sii impair et 7 pair, 
h=UL2sn 
27 LÉ 27 sii impair et 7 pair, 
| cos’/(0)sin/*?(0) cos((2h + 1)0)d0 = HO nn 
ÿ 187 sii impair et 7 pair, 
LR DORA 
0 si non, 
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[ cos’ (0) sin/*?(0) sin((2h + 1)0)d0 — 


[ | cos’ (8) sin/*?(9) sin((2h)0)d0 = 


le cos’ (0) sin/*1(9) sin((2h + 1)0)d0 — 


[ cos’(0)sin/*1(0) cos((2h + 1)0)d0 = 


| : cos’(0)sin/*1(0) sin((2h)0)d0 — 


[0 sin/*1(0) cos((2h)0)d0 = | 
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TA sit pair et 7 impair, 
A STE 

API sit pair et j impair, 
RE 

rA?# + sit pair et j impair, 
h= 012.4 

0 si non, 


nC2 sii impair et j impair, 


h= 02: 
rC2? sii impair et 7j impair, 
RD Lou 
0 si non, 


mD'+1 sit pair et j pair, 


ho 0. 192250 
TD sit pair et j pair, 
= 0,1:22: mm 
mDES? sit pairet j pair, 
DOS 
0 si non, 
RE sit impair et j impair, 
be? ssn 
RE si impair et j impair, 
RQ en 
mÉPRTI si impair et 7 impair, 
he DL Dot 
0 si non, 


mF?*  sii pairet j impair, 


A OP Te 
r2r sit pair et j impair, 
Ds DL. 
ah sit pair et j impair, 
16 90 ON D 
0 si non, 


nG?* sit impairet j pair, 


R= UE 220 
rG?i sit impair et 7 pair, 
6 D PU EP 
0 si non, 


TA sii impair et ? pair, 
hi= is 2:20 
07 moe sit impair et j pair, 
} cos *1(0) sin/*1(9) sin((2h + 1)0)d0 — he 0,162; 0m 
ô TH sii impair et 7 pair, 
BU, al 
0 si non, 


2h+1 ÿ à : 2 3 ; 
rl; sii pair et j impair, 


AVR 
2 LE Fi sii pair et j impair, 
| cos*1(0) sin/*1(9) cos((2h + 1)0)d0 = OT sn 
ù La Cd ifà pair et j impair, 
LRQ Pi PL 
0 si non, 
TKE? sii impair et ? impair, 
ROLE 
2x mA sit impair et j impair, 
| cos" *1(0) sin/*! (0) cos((2h)0)d0 = = 0"1,2%::.m 
ù rK2 sii impair et ? impair, 
À ON 
0 si non, 


2h CS . . . 
ml; sit pairet j pair, 


M | , 0 A EP EEE 
J cos" T"(4) sin/**(9)sin((2h)0)d0 = 4 mL sii pairetj pair, 
0 ot 

0 si non. 


2h+1 A2h+1 AÂ2h+1 p2h+1 H2h+1 B2h+H1 2h 2h p2h+1 
Ou Qij Bis dés Ai; , A;; À; , B;; , B;; , B;; Ce C7 pe ? 
L h+1 fr2h+1 
J Hi 


Y 


y 2h41 2h41 F2h4I F2h4+1 pr2h 2h fe2h 2h T 2h 
HET IE LS TE, Kf7, Kiÿ, Ki, Li et Li sont des constantes non nulles. 


AMijxn  Sii impair et j pair, 
k pair et h impair, 


27 
| cos TE+1(0) sin/##3(9)d0 = 4 rMijn sit pair et j impair, 
8 k impair et À pair, 


0 si non, 
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ANijkh sii impair et ? pair, 
k est pair et h est pair, 


oui pair et 7 pair, 
k impair et À pair, 


27 
| cos! +1(9) sin/}42(0)d0 = 4 TNijun Sit impairet j impair, 
0 


k est pair et h est impair, 


ANijkh sit pair et j impair, 
k impair et h impair, 


0 si non, 


REte sii impair et ? pair, 
k impair et h pair, 


T À; jun sii pair et 7 impair, 
27 : é à 
À cos! #(0) sin/t#+4(0)d0 = PAAER EMRDAE 
0 


APijxh  Sit impair et j impair, 
k impair et À impair, 


0 si non, 


TQijkh sit pair et j pair, 
k impair et À impair, 


ou? impair et 7 impair, 
k pair et h pair, 


27 
1 cos +#+1(0) sini+#+1(0)d0 = À rOQjun sii pair et j impair, 
ê k impair et h pair, 


ou impair et 7} pair, 
k pair et h impair, 


0 si non, 
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27 
| cos! T# (0) sin) ##3(0)d0 = 
0 


27 
J cos! #29) sin }+2(9)d0 — 
0 


27 
| cosiT##2(9) sinft*+1(9)d0 — 
0 


A Rijkh sii pair et ] pair, 
k pair et h impair, 


ou? pair et ? impair, 
k pair et À pair, 


TEE sii impair et 7 pair, 
k impair et À impair, 


or à impair et j impair, 
k impair et À pair, 


0 si non, 
Tien sii pair aet } impair, 
k pair et h impair, 


TTijen sit impair et j pair, 
k est impair et h est pair, 


TER sii impair et ? impair, 
k impair et À impair, 


0 si non, 
AÜijkh sii pair et } pair, 
k pair et h impair, 


ou? pair et ? impair, 
k pair et h epair, 


RÜijrkh sit impair et j pair, 
k impair et h impair, 


ou à impair et 7 impair, 
k impair et h pair, 


0 si non, 
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A Vijkh sii pair et} pair, 
k impair et À impair, 
7 Vjrn sii pair et 7 impair, 
27 : : : 
;. cosi+*#1(9) sini+*+3(9)d0 — k impair et h pair, 
0 


ou? impair et ] pair, 
k pair et h impair, 


0 si non. 


Mine, Nes, Ni Noirs Qui, Qins Potns Pins Pérr, Roses R17KR, T4, 


Lib bpevens Ve VE 


sont des constantes non nulles. 
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